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Abstract. The geometric theory of vortex tunnelling in superfluid liquids is developed. 
Geometry rules the tunnelling process in the approximation of an incompressible superfluid, 
which yields the identity of phase and configuration space in the vortex collective co-ordinate. 
To exemplify the implications of this approach to tunnelling, we solve explicitly for the 
two-dimensional motion of a point vortex in the presence of an ellipse, showing that the 
hydrodynamic collective co-ordinate description limits the constant energy paths aUowed for 
the vortex in configuration space. We outline the experimental procedure used in helium II to 
observe tunnelling events, and compare the conclusions we draw to the experimental results 
obtained so far. Tunnelling in Fermi superfluids is discussed, where it is assumed that the 
low energy quasiparticle excitations localised in the vortex core govern the vortex dynamical 
equations. The tunnelling process can be dominated by Hall or dissipative terms, respectively 
be under the influence of both, with a possible realization of this last intermediate case in 
unconventional, high-temperature superconductors. 
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1 Introductory considerations 

The primary notion we have of a superfluid is that it shows no dissipation of the 
flow, under certain, well-defined conditions. Friction can be caused by the creation 
of elementary quasiparticle excitations in the superfluid, that is, irreversible energy 
transfer from the coherently moving superfluid state to incoherent degrees of freedom. 
This can occur if the superflow (relative to some reference frame) exceeds the Landau 
critical velocity, creating the excitation, and thereby reducing the superfluid current. 
A different kind of dissipation can be caused by a topological excitation, the quan- 
tized vortex, representing a travelling defect structure in the order parameter of the 
superfluid. The dissipation mechanism is then represented by the vortex crossing the 
streamlines of the flow, diminishing the superfluid current by reducing the superfluid 
phase difference between points on a line perpendicular to the vortex motion. To 
cause current reduction, a vortex first has to be generated. The task of this paper 
is to develop a formalism describing a quantized vortex entering a superfiuid at the 
absolute zero of temperature by the quantum mechanism of tunnelling. 

The nucleation theory of quantized vortices in the Bose superfiuid helium II has 
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been an elusive subject ever since the existence of quantized vortices was conjectured 
by Lars Onsager in 1949 §. The difficulty to grasp their coming into existence in 
a quantitative manner from first principles has one fundamental reason: There is 
no microscopic theory of this dense superfluid. We do not know how to describe the 
motion of a vortex on atomic scales, where this motion is governed by the full quantum 
many-body structure of the superfluid. One can even go further, and state that we 
even do not know precisely what a vortex should be on these small scales. Its very 
definability as a stable topological object essentially depends on the usage of a large 
scale approach. 

Before we describe the problems inherent in vortex nucleation theory, to clarify 
terminology we use here and further in this work, it is advisable to fix some notions. 
The treatment of a fluid will be called hydrodynamic or large scale, with the frequency 
of excitations of the fluid approaching zero for large wavelengths, if the underlying 
atomic structure of the fluid is not relevant for the phenomena under investigation. The 
fluid is reasonably well, on these large scales, approximated to consist of structureless, 
pointlike particles. In particular, a flow field can be deflned, with the prescription that 
a volume element moving with a certain flow velocity contains enough particles for the 
hydrodynamic formulation to make sense. 

The definition of semiclassicality is related to the existence of quantum 'fluctua- 
tions' or, better, the fact that a quantum mechanical variable is indeterminate in its 
value with respect to the outcome of different measurements. If quantum fluctuations 
are small against the expectation value of a quantum operator, we speak of the semi- 
classical limit. The notions of hydrodynamic and semiclassical in the dense superfluid 
helium II are mutually corresponding to each other, and a hydrodynamic treatment 
on macroscopic scales has semiclassical accuracy. A mean field theory asserts that 
it is reasonable to replace the fleld operator by its expectation value, in particular 
within expressions of order higher than the second in the operator and its conjugate 
(an example is to replace two of the field operators in the quartic interaction term of 
the second quantized Hamiltonian by their expectation value). It therefore makes the 
assumption that the system can be treated on any scale by making use of such a re- 
placement. This procedure certainly cannot be used in the case of helium II. Though, 
on large scales, this dense superfluid is, to a good approximation, described by an 
order parameter function (j), on the atomic scales of order the coherence length ^, it 
is necessary to solve the full second quantized problem. There the indeterminacy of, 
e.g., the operator particle density is as large as its expectation value. 

A measure of the applicability of mean field theory can be obtained if we consider 
the effective strength g of interaction. It can, in the limit of long wavelengths, be 
defined to be the spatial integral of the two-body interaction of "'He atoms, and is 
related to the s-wave scattering length a the speed of sound Cs, and the bulk 
density po hy g = ATrfi^a/m = mc^/po (repulsive interaction, a > 0, [^). A mean 
field treatment is useful if poa'^ ^ 1 which implies a = Tr{mcs/h)'^ / po <C d, with 
d = Pq the interparticle spacing. Using the data of |^], we have in helium II (at p ~ 1 
bar), d ~ 3.58 A, Cs — 238 m/s and thus a ~ 8 A, increasing further to approximately 
double this value at solidification pressure. The scattering length, that is, the effective 
range of the scattering potential of the atoms for long wavelength excitations, is larger 
than the interparticle spacing, and the condition for the applicability of mean fleld 
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theory is violated. This fact entails that one has to consider helium II as a strongly 
coupled system, which has caused quite formidable computational efforts trying to 
understand its behaviour on sub-a scales Based on the above considerations, 

we will understand a dense, strongly coupled superfluid as one fulfilling the condition 
a > d and, correspondingly, a dilute, weakly coupled superfluid as one having a <C 

d, the Bose-Einstein condensed atomic vapours being examples of this kind of 
superfluid. The condition a > d implies that, in strongly coupled superfluids, the 
coherence length is of order the interparticle spacing. 

We have already mentioned at the beginning that there are two classes of exci- 
tations related to dissipation, the first being Landau quasiparticles, the second one 
defects in the order parameter. There arises the question if these different branches 
of excitations merge in a certain manner, approximately for a wavevector 2TT/d cor- 
responding to the interparticle spacing. The roton, i.e. the excitation at the local 
minimum of the excitation spectrum of superfiuid helium II, has a wavevector kr 
which is very close to this value, kr ~ 27r/d. The idea that the roton could correspond 
to a vortex ring of the smallest possible size (namely, such that just one atom can pass 
through it), goes back to the seminal papers of Feynman §-[|lO|. Evidence for such 
an identification would indeed provide justification to speak of some kind of, at least, 
similarity on the scale of the interparticle distance, of the two types of excitations. 
There is still a debate going on about this possibility (recent contributions are found, 

e. g., in ||3^,]3^]) In this work, we will, however, not enter into such a discussion, for it 
is evidently still a long way to a complete understanding of the microscopic dynamical 
behaviour of vortices, as even their stationary microscopic character (provided that a 
vortex ring of atomic size is properly definable at all), is not yet completely clear. We 
will take the (classic) point of view that a vortex is a topological object, well-defined 
as a defect structure in the order parameter. Then, a vortex ring of atomic size is 
termed virtual, because it does not constitute a topological object on this scale. 

The lack of a microscopic idea of vortex motion makes it necessary to resort to a 
hydrodynamic theory of the motion of a vortex, be this motion in real or imaginary 
time, the latter of interest for the tunnelling processes we intend to investigate. In a 
hydrodynamic treatment, the existence of the vortex as a semiclassical object has to 
be assumed ab initio. No details of the underlying microscopic dynamics, i.e. of the 
actual nucleation event, are to be described in such a theory, but only the laws which 
rule vortex motion on curvature scales well beyond the atomic one. The microscopic 
dynamical behaviour of the vortex is, in such a description, bound to appear only 
in cutoff parameters determining the borderline to the microscopic realm. We will 
consequently consider the vortex motion as it results from the Lagrangian in terms of a 
collective co-ordinate for the vortex, which is useful as long as the curvature radii of the 
line described by this co-ordinate are much larger than ^. Such an approach in terms 
of a collective co-ordinate makes sense and is physically meaningful, if we additionally 
assume that the potential barriers, through which the vortex has to tunnel, have 
themselves effective curvature radii well beyond a scale O(^). If tunnelling events can 
be actually observable, then, depends for a given superfluid at a given temperature 
T on the ratios £,/d, T*/T, where ksT* = h^/{md'^) is a characteristic quantum 
energy of the quantum fluid constituted by particles of mass m and T* an associated 
temperature; for helium II, T* w 1 K. This dependence can qualitatively be understood 
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as follows |ri[ . Consider a vortex ring of radius Rq, which has in the bulk an energy 
E{Rq) = {mpoK^ /2)Ro\n{Ro/ Rc) (see section 2.5), where Rc — 0(0 is the vortex core 
radius. The relevant quantity to compare this energy with, is the thermal energy ksT. 
If E{Ro) 3> ksT, quantum tunnelling is exponentially suppressed. Writing the energy 
in terms of the above quantities, E{Ro)/kBT — {T* /T){Ro/d)lii{Ro/ Rc) (barring a 
factor of 1/4 for pair-correlated superfluids). The ratio (,/d then effectively enters, 
because the smallest possible value of the radius of the vortex appearing in the fluid 
is Rq > Rc = 0{£,). Considering the fact that thus simultaneously T ~ 0(T*) and 
Rq ^ 0{d) are to be fulfilled, helium II is the most promising candidate for quantum 
tunnelling to happen. The conventional superconductors and '^le, with their large 
T*/T as well as £,/d are, already on this ground, ruled out. Hence the only possible 
candidates remaining for an observation of quantum tunnelling of vortices are, save 
for helium II, high-Tc superconductors. 



Overview 



To give the reader a concise impression of what follows, we provide here an overview of 
the principal directions to be pursued, and ideas to be developed in the three sections 
of this work to follow. The theory of quantum tunnelling is developed in the next 
section. It is shown that in the limit of long wavelengths, which we are required to be 
using in a dense superfluid, the probability of quantum tunnelling is predominantly 
determined by external geometrical quantities connected to the geometry of the flow. 
The tunnelling exponent is separated into a dominant volume contribution solely as- 
sociated with the tunnelling path of the vortex, which is a contribution independent 
from the fact that the fluid is compressible or incompressible, and a subdominant area 
contribution associated with that same path anc? the fact that the fluid is compressible. 
The dependence of the dominant volume term in the tunnelling exponent on geomet- 
rical quantities is exemplified by the analytically solvable problem of a point vortex in 
the presence of an ellipse, where the long wavelength, collective co-ordinate treatment 
is shown to impose strict constraints on the motions allowed for the vortex. In the third 
section the experimental procedure to observe the temperature independent quantum 
mechanical triggering of vortex generation, which we wish to explain, is demonstrated. 
We discuss the data obtained in these experiments, with particular emphasis on the 
applicability of our predictions in this work. The fourth section describes some as- 
pects of pair-correlated, charged Fermi superfluids. We discuss, in particular, the role 
which might be played by the existence of bound quasi-particle states for observable 
vortex tunnelling phenomena in superconductors. The high-Tc superconductors, on 
account of their small coherence lengths, play a prominent role in these considerations, 
as emphasised above. It is explained that, even for very low temperatures, in the case 
of unconventional (d-wave) high-Tc superconductors of practically feasible purity, the 
tunnelling phenomenon is not adequately described by the theory of sections one and 
two, as this were the case for conventional (s-wave), extreme type II superconductors, 
on a length scale well below that of the magnetic penetration depth. 
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2 Quantum tunnelling of vortices 

2.1 General introduction 

The quantum mechanical phenomenon of tunneUing has attracted attention from re- 
searchers in theoretical and experimental condensed matter physics, field theory and 
other areas. It belongs to the most remarkable properties solely pertaining to systems 
obeying the laws of quantum mechanics. Essentially, quantum mechanical tunnelling 
in a condensed matter system is described by the motion of a few degrees of freedom 
(subsumed in what follows into the term 'particle') under a potential barrier in con- 
figuration space. The particle is assumed to have less energy than represented by the 
height of the barrier (setting the bottom of the potential equal to the zero of energy). 
Because of wave-particle duality, a wave function can be associated with the particle. 
This wave function is damped under the barrier or, in other words, the particle trav- 
els with imaginary momentum there. To some extent and with some probability, the 
particle is thus located under the barrier. It can then even completely penetrate it, if 
the damping is small enough, getting with some nonzero probability from one side of 
the barrier to the other, whereas classically this is a completely forbidden process. 

There are different means to describe the tunnelling motion. The canonical way 
is to calculate wave functions for the problem, solving, e.g., the Schrodinger equation 
for the potential of interest with appropriate boundary conditions imposed |p^ , |l3| . 
The most popular formal means to investigate tunnelling, however, is provided by 
the calculation of the Euclidean action of the system along the tunnelling trajectory 
[H, |l5|. The Euclidean action is obtained from the Minkowski action by rotating 
to purely imaginary times Se = ~iS[t — > —it,,]. The time on the imaginary axis of the 
complex i-plane will be denoted te. If the action is dominated by the classical path 
in Euclidean time, the tunnelling probability can be calculated in the semiclassical 
approximation. The corresponding solution of the second order Euclidean classical 
(field) equations of motion with finite action is called instanton. The name stems 
from the fact that the instanton is a particle-like object localised in Euclidean time. 
It exists, so to speak, just long enough for the actual particle to tunnel. 

In the semiclassical limit, the tunnelling probability for a given energy E is taking 
the form 



In this relation, Se{E) is the Legendre transform of the Euclidean action Sc{Te) ^ h 
as a function of Tg, the Euclidean period of motion H: 



P{E) = A{E) exp 



Se{E)^ 



(1) 



h 



Se{E) = SeiT,) ~ — f Te = Se{Te) " ETe . 



(2) 



The quantity A{E), frequently called pref actor, represents the influence of fluctuations 
around the classical path. It is essentially given by the inverse determinant of the co- 
efficients of the second order deviations from the classical path in the action IgI EbI . 
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2.2 The collective co-ordinate action of the vortex 

It is a quite commonly accepted wisdom that any complex condensed matter prob- 
lem remains intractable if we do not single out certain central, collective degrees of 
freedom, termed in general 'collective co-ordinates'. This is possible because there 
are conservation laws and symmetries governing the behaviour of the system as a 
whole: We can actually describe essential features while not referring explicitly to the 
10^" • • • 10^^ particles and their interaction. 

The obvious choice for the vortex collective co-ordinate is its center X'^{t, a), which 
also indicates the center of topological stability and thus conserved topological charge. 
That this co-ordinate represents the vortex sufficiently accurate requires that we con- 
sider scales much larger than the vortex core size of order the coherence length ^. 
Furthermore, we assume that there is a canonical collective vortex momentum related 
to this central co-ordinate. The action (||) then is 



SeiE) = -i I j)dtedaX ■ P = -i (jxlaj) dX ■ P 

dK P, (3) 



where we defined the imaginary differential co-ordinate vector dK = —idX of the 
vortex.^ The parameter a labels points on the vortex string, and P is the canonical 
momentum per cr-length. 

Note that the co-ordinate differential vectors in (|^) are no function of cr, as the 
co-ordinate position vectors themselves of course are. The closed time integral indi- 
cates that we take the integral over a full period of the motion. That such a periodic 
motion exists, is of course a highly nontrivial assumption for arbitrary dimension of 
the phase space. Only in an effectively one-dimensional problem (one spatial dimen- 
sion), respectively for multidimensional systems separable into such one-dimensional 
problems (cf. |^ §48), closed phase space trajectories necessarily exist. 

2.2.1 Contributions in the tunnelling action 

In two dimensions and in a conventional, electrically uncharged superfluid, vortices 
and charged particles have identical dynamical equations in the hydrodynamic limit. 
According to the three-dimensional extension of this duality, i.e. as a vortices-charged 
strings analogy, the vortex Hamiltonian takes the form |0] 



Hv = 



j)d(7^ 



M,cl^-\-{P^qaf+'-^Q' 



da ( — 
2 



(4) 



The first integral represents the self energy of the vortex, which has static, kinetic, and 
elastic contributions, respectively. The arc length of the line is written as y/^da, and 
the vector Q is a perturbation of some equilibrium string configuration perpendicular 
to the tangent vector X' of the line. The rest frame mass Mq is given by Mq — Eq/cI 

^K" is not to be confused with a wave vector. We could have chosen as well to incorporate the —i 
into the (imaginary) momentum P. Crucial is only that Sc is a real quantity. 
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[ p8[ , where Eq is the logarithmically divergent static self energy of the vortex (per 
unit length). In helium II this energy reads Eq = [N^ K^mpo / {An)] In (8Rc/£,e'^) , with 
po the bulk number density, m the helium mass, k = h/m the velocity circulation 
quantum, Rc the infrared cutoff (the local curvature radius of the line) , ^ the ultraviolet 
cutoff (the core size), and C a constant characterising the core structure. The second 
integral represents the 'electrostatic' energy, in which we split the scalar potential 
a" into a Coulomb contribution from the interaction with other vortices, and the 
interaction with a nonvortical background flow, a^. The vectorial generalization of 



the stream function ip of classical hydrodynamics |27 , is related to ao = —a" by 
ao = poip ■ X' . The factor 1/2 in front of stems from the Coulomb gauge for the 
vector potential, diva = 0, and makes sure that the energy of each vortex is counted 
only once. The 'charge' q = Nyh is given by the winding number Ny multiplying 
Planck's quantum of action.^ 

The gauge potentials are derived from the external flow field at the position of the 
line element at cr by a gauge invariant duality relation fl^ . In its nonrelativistic form 
needed here, this reads 

rota = -pX' , (5) 
dta + Va" ^ X' xj , (6) 

where j = pv is the conserved particle current, dtp + divj = 0. The quantity of crucial 
importance in the Euclidean action of constant energy (^) , determining the tunnelling 
exponent, is the canonical momentum 

p ^ pine ^ pkin ^ ^ Mo^X . (7) 

It consists of a contribution P™'^, which is related to the Magnus force acting on the 
vortex, and a second contribution P*^'" related to the existence of a nonzero vortex 
mass, that is, to a finite compressibility and thus finite Cg- It follows by integration of 
(^ that the ratio of the momentum contributions i-**^'", P'"'^ is in order of magnitude 
~ (NyK) / {cs\X\){\X\/cs) (neglecting the vortex energy logarithm). Hence for large 
scales (large curvature radii), and small velocities, that is in the hydrodynamic limit, 
P""'" is dominated by P'"^ 

The Euclidean action is given by 

+q J (j) dtedcr - ia ■ X . (8) 

In case that a tunnelling process of constant energy is under consideration, the quantity 
of interest is the action as a function of constant energy (||). According to the relation 
(^, this action consists of a part related to the vector potential and another part 

^ In pair-correlated Fermi superfluids, q = Nyh/2, where the number density in and ^ is 
understood to refer to the "elementary" particles constituting the superfluid, and not to the Cooper 
pairs. 
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related to the vortex effective mass {X = Q) 

dtpdaX ■ 





pine pkin 


j <j) dtedaX ■ 







-iqa + Mo^X\ . (9) 

The EucHdean action sphts into a part S^J^"^, due to the interaction of the vortex with 
an (approximately) incompressible background superfluid, and a part 5*^™ which can 
be ascribed to the kinetic ('vortex matter') term in the vortex momentum. We will 
now show that S"^'^{E) is given by a volume associated with the path the vortex line 
traces out in configuration space, whereas S^^'^ is connected with an area associated 
with that path. 

To demonstrate this, we integrate relation (^, multiplied with q, to obtain (p — po) 



P'^'^dX^da = NyhpQ j j j y^dX'dX-^da , (fO) 

where g is the determinant of the coordinate basis on the line (unity for a triad) . The 
closed surface with surface elements of magnitude dX'^da {a— 1,2 is the index of the 
two Q-directions), encloses the total volume with local element ^ dX^ dX^ da traced 
out by the line on its path. Further, using the gauge freedom for the momentum, we 
can express the action S™'^{E) by the volume integral 

SfiE) ^ 2^^^ po III y^dZ^dZ^ da , (11) 



wherein the co-ordinate differentials are defined to be 

dZ^ = cosa dK^ + sina dK^ = —i (cosadX ^ + sinadX^) , 

dZ^ = -sinadX^ +cosadX^ . (12) 

The angle a{a) parameterizes in these differentials the local (co-ordinate) gauge free- 
dom for the momentum, of rotations about the line tangent X' . It expresses the fact 
that one degree of freedom is still available, namely that for the direction of the local 
gauge dependent momentum P^^'^{a), even after a local basis on the string has been 
chosen. For the components of f '"'^ in the two Q-directions Cj and eg the relation 

d^P'r-diPr^N^hpo^ (13) 

obtains |26 , cf. relation (|^). The simplest example is a rectilinear line in z-direction, 
for which the local momentum can rotate in the x-y plane. The gauge invariant quan- 
tity is the integral §§ P™^dX°'da, which is left unchanged by the rotation freedom. 
The relation (|l||) implies that in the hydrodynamic Hmit of |P'''"|/|P'"''| 0, phase 
space and configuration space become indistinguishable, since the momentum compo- 
nents then become functions of the co-ordinates alone, and are no longer independent 
variables | p^ . 

It is illuminating to go back to our 'electrodynamic' quantities and rewrite 

S'^"'{E)^iq 1 1 1 B^^dX^dX^ da , (14) 
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where Bo- = — po is the (nonrelativistic) 'magnetic' field, pointing antiparallel to the 
direction of the hne tangent. The part S™''{E) is thus the Aharonov-Bohm type Berry 
phase [ p^ in the Euclidean wave function of the adiabatically moving quantum object 
vortex. 

The part of the action S^"^ explicitly involves the vortex dynamics, and thus can 
not be calculated by knowledge of the vortex path geometry alone, as this was possible 
for the part S]^'^. Treating the influence of the mass as a small perturbation on the 
vortex path, the ratio of the actions is of the same order. More exactly. 



in 



H---)iQa/c,)^dX^da 



- '^Po^^^^^ ffiQa/cs)VldX''da (A^, = 1, He II). (15) 

The last line is valid for a unit circulation vortex in helium II, in which the approximate 
equality 

K ~ Cs2^ (in helium II) (16) 

holdsj^ The logarithm is assumed to vary only slightly during the tunnelling process, 
so that it is written in front of the integral. In the length scale domain of interest, 
ln(- • ■)/2tt = 0{1). The dots, • • •, indicate an average over the argument of the vortex 
energy logarithm. 

Essentially, ( p^ tells us that the kinetic part of the Euclidean action in units of h 
depends on an area (element) multiplied by the (local) velocity of the vortex in units 
of the speed of sound. The coherence length ^ times this area gives a volume, which 
multiplied by the bulk number density finally yields a dimensionless action. The full 
Euclidean action (^) thus takes the schematic form (for simplicity again displayed in 
the unit circulation case of helium II) 



SeiE) 



^ Po(2vrl]('^)+ln(...)eS('')[aX/(c,9i)]) . (17) 

The volume rj^'*) = -i fff ^/g dX^ dX^da. The effective surface 

S^'^) = j^j^{Qa/cs)VldX''da (18) 

obtained by integrating over the surface enclosing Q^'^'^ is a functional of the vortex 
velocity scaled by the speed of sound (d indicates the spatial dimension) . 

To summarize, in a dense, strongly coupled superfluid, the first term dominates 
the second term in ( p7[ ) for the following reasons: 

'^This relation assumes that ctlj < I; < a, with ^ nearer to the lower bound (a is the interparticle 
spacing, ctlj = 2.556 A the Lennard-Jones parameter of the atomic helium interaction), which is 
consistent with quantum many-body and density-functional calculations [ p^ , |34| ^ . 
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a. The large scale, collective co-ordinate limit requires that the scales to be consid- 
ered are much larger than ^. The corresponding volumina and areas have to be 
very much larger than ^'^ and respectively, and hence fi'-''' ^ S^T,^'^^ . 

b. The area contribution of S'g'" is additionally suppressed by the vortex velocity 
divided by the speed of sound, i.e. T,^'^^ ^ dil,^'^\ where cJil^''^ is the proper 
surface area enclosing fl^'^K 

This dominance of S'^"'^ over S^™ is in contrast to the case of relativistic (string) objects 
moving with speeds of order Cs Hi], |2^. Under this circumstance, S'^'" is of the same 
order as S™^ and the action is of order Se/fi « po ln(- ■ ■) S, dn'^'^\ If one neglects the 
dependence of the logarithm on the co-ordinates, this is essentially the Nambu action 
((2^), in units of h, up to a factor of order unity. 

We assumed that the vortex path in phase space is closed. As a consequence, the 
number of particles in the effective volume on the right-hand side of ( [T7| ) is quantized 
according to 

Se{E) = {N^'^^ + a)h 4=» S^{E)/h = 2TT{M'''> +a). (19) 

The number of particles in the effective volume (including the small kinetic contribu- 
tion on the right-hand side of (p^)), plays the part of the Bohr-Sommerfeld quantum 
number in semiclassical quantization. The number a is of the order one and signifies 
the onset of the microscopic quantum regime. In the semiclassical approximation, 
N'-'^'> ^ a must hold, so that A^^**) ^ 1 gives, as usual, a direct measure of semiclassi- 
cality. 

2.3 Geometry of Quantum Tunnelling 

2.4 Galilean in variance violation 

At the absolute zero of temperature, a homogeneous nonrclativistic supcrfluid has 
Galilean invariance, that is, physical contents are invariant under co-ordinate trans- 
formations to any reference frame moving at constant velocity. If we approach absolute 
zero, which is what is actually realized in experiment, we expect the tunnelling rate 
to make no abrupt change as the temperature is lowered. Thus the result for the rate 
we obtain at T = should also be valid for temperatures slightly above zero (we will 
estimate the temperatures, for which this is no longer the case, later on). 

Because we can always transform to the rest frame comoving with the superfluid, 
the tunnelling probability at T = equals zero if Galilean invariance remains unbroken: 
In the rest frame there is a tunnelling barrier of infinite height, the logarithmically 
diverging vortex self energy. Hence it is necessary to explicitly include the violation 
of Galilean invariance by a flow obstacle into any calculation of tunnelling rates for 
Galilei invariant superfluids at absolute zero. The necessity of invariance violation for 
tunnelling to be energetically allowed thus stems in the real superfluid from the fact 
that it is possible to invariantly transform to the rest frame of the superfluid, for any 
allowed velocity of flow. 

Having thus shown that it is strictly required that Galilean invariance be violated, 
we are in demand of constructing an explicit solution of the following problem. Given 
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a vortex in the presence of an invariance breaking flow obstacle, we have to calculate 
the Hamilonian energy Hy of the vortex in the Hamiltonian as a function of the co- 
ordinates, by solving one of the equations of motion corresponding to this Hamiltonian 
respectively the Euclidean action (^). We then set Hy equal to a constant E, to 
calculate the vortex trajectory of constant energy, which finally yields the action (||) 
respectively (^), in the form of (17). It is obvious that without a very high degree of 
symmetry, this is a task necessitating a quite formidable computational effort. 

We have seen that the problem of determining the dominant contribution S^"'^ is a 
geometrical problem, because the phase space co-ordinates are functionally dependent 
on the configuration space co-ordinates. We will therefore restrict the discussion in 
what follows to cases which elucidate in particular the geometric nature of quantum 
tunnelling, that is, concentrate on the behaviour of S'^"'^, which is a functional of 
geometrical quantities. 



2. 5 The vortex half-ring case 



To flesh out the discussion which has been so far quite abstract, we now give some 
quantities relevant in the tunnelling problem for a vortex. Imagine, for concreteness, 
a singly quantized vortex half-ring of radius R with its circulation axis standing 
perpendicular on the plane y = 0. Over the plane is flowing liquid with a velocity u at 
infinity from right to left (in the negative x-direction) . Consider, first, the case that 
the invariance-breaking asperities on the surface are small-scale (it will become clear 
in a moment what 'small-scale' means). Neglecting the elastic and kinetic terms in 
the Hamiltonian, as well as perturbations of the fiow generated by the asperities, the 
conserved vortex energy is simply 



Hv = hpQ / d(j) 



i? In 



An 

If we normalize the energy by 
Anm 



8R 



1 



and solve for the path of constant total energy Eq, we get the relation 

8i?/C 



R 



2nu 



In 



exp{C + Eo/{nRj) 



(20) 



(21) 



(22) 



This equation has two solutions for the radius R. One of them is, in the case of 
u <C K/2n(, and small E^/nR ^ C, located far away from the surface cc = 0, at 



Rq 



2nu 



In 



8K/(27rO\ 
exp(C) ) 



(23) 



The other solution is a half-ring of radius in the order of the coherence length, R^ — 
O(^) (which will not be of further interest here; it signifies the path of the (virtual) 
vortex trapped at the flat boundary). 

The closed path in phase space, needed to evaluate (^), can be obtained as follows. 
Assume that the small asperities are approximately of the shape of oblate rotation 
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ellipsoids with half-axes b perpendicular and a parallel to the flow and define the 
sharpness (3 = b/a > 1. Then, the half-ring trajectory over such a small ellipsoid, 
coUinear with the ring axis, were given by 

= b'^ + 13^ K'^ , (24) 

using the complex co-ordinate K = —iZ, and provided that the vortex exactly follows 
the surface (up to a constant 0(^)). This trajectory hits the one far away from the 
surface, given by Rq =const., at Kq ~ Ro/P (provided that b <^ Rq). The closed phase 
space trajectory thus begins at the ellipsoid top {K = 0, i? ~ 6), then propagates along 
the line given by (p4[), meets the constant Rq at Kq, then follows this constant to —Kq, 
and finally follows the branch of (|2^) for negative K back to the ellipsoid top. 

The tunnelling exponent, using the gauge Px — (l/2)/i/5o^^ for the momentum, is 
thereby 



S'riEo) = 2hpo f dcf, f "dKliRl-P^K^) 
Jo Jo ^ 



= hpoYRo/P^hpon<^^K (25) 

The tunnelling volume and thus the exponent may consequently be reduced and the 
tunnelling probability enhanced if a surface with sufficiently sharp peaks perpendicular 
to the flow is present. In particular, if we were able to reduce the small half- axis 
to a > ^, and still have 6 a, we could reduce the tunnelling volume to fl'^'^^ ~ 
(27r/3)i?Q ^(i?o/^)- It is worthwhile to point out that this reduction is not due to an 
enhancement of the flow velocity at the ellipsoid top: For any value of /3, because of 
cylindrical symmetry, the velocity at the top is exactly 2u, just as for the half-sphere. 

The trajectory ( p^ is correct only in lowest order of /3 — 1. In particular, for /? = 1 
(a = b), the half-sphere, the result for the tunnelling exponent ( p5| ) is exact in the 
low velocity limit. This can be shown by solving the problem for the sphere exactly, 
which has been done in |23] . The result for the stream function part of the potential 
in the Hamiltonian (^), a" = —poip, corresponding to Stokes' stream function -tp [ p7| , 
has been checked by the author, and is in Coulomb gauge 

l^c + ^.=R^^Q.M^)-luR^[i-{-^^)') , (26) 

with ui = 1 + (1/2)[(Z^ + _R^ — a^)/(ai?)]^. The function Q 1/2 is a Legendre function of 
the second kind (|^^, No. 8.821), a solution of the Poisson equation At/i = kS^^^ {x—X) 
for the combined spherical-cylindrical symmetry of the problem at hand. 

The problem of a undeformed, massless half-ring situated at a half-sphere, with 
their axes coinciding, is actually the only nontrivial problem in three spatial dimensions 
solvable with reasonable effort analytically. It has just the maximal symmetry in three 
dimensions. It does not, however, contain the crucial ingredient for a full geometrical 
analysis of quantum tunnelling: The curvature of the flow obstacle is constant. The 
result that, for the sphere, the vortex half-ring follows exactly the surface (up to a 
constant distance of order ^), has led in psf to the assumption that it is admissible 
to continue this result to the ellipsoid case, in the manner of (E4). We will show 
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in the subsection to follow that this is incorrect: The analytic continuation of the 
sphere result to the ellipsoid case of arbitrary (3 is incompatible with the geometric 
requirements imposed by a hydrodynamic, collective co-ordinate treatment. 

2.6 Analytical solution in two dimensions 

A treatment of our (boundary) problem in two spatial dimensions is advantageous 
for the following reasons: 

i. The additional spacelike co-ordinate, namely the arc length parameter, compli- 
cates the analysis enormously because of the locality of any variable in a and 
the existence of elastic energy. 

a. Even if we treat an undeformed ring or half-ring, i.e. neglect elastic energy 
and locality altogether, any 3d problem which has not the maximal symmetry 
described above, requires a multitude of image vortices and gets quite intractable. 

in. In two dimensions, we have available the tools of conformal transformation, 
allowing for comparatively simple calculations of analytical solutions. 

iv. The geometrical features, which are dominant in the hydrodynamic, collective 
co-ordinate limit, as expounded at length above, are most clearly seen in two 
dimensions. 

2.6.1 The solution for the circle 

The basic solution from which we start is that for a vortex in the presence of a half- 
circle of radius d at an otherwise flat boundary (cf. Fig. |l]). The complex plane of this 
original solution is called the Z-plane (the uppercase letter does here not imply that 
a vortex position is meant). The imaginary part of the complex potential [ p7[ gives 
the stream function tp = ^[w], whereas the real part is the usual velocity potential. It 
follows that a single vortex at Zi has complex potential w{Z) = — i(K/27r) ln[Z — Zi]. 

The boundary condition to be fulfilled is obviously that there be no flow into 
the surface consisting of the line Y = Q and the half-circle. This amounts to the 
requirement that the line F = (for |y| > d), F = Vd^ - X"^ (for |y| < d) is a 
streamline of constant ip = 0. Such a requirement can be met by using the technique 
of image charges quite familiar from electrostatics: Our vortex problem is completely 
equivalent to that for a 'charge' situated near a perfectly conducting surface, with no 
tangential 'electric' field. 

The complex potential generated by the image vortices and acting on the vortex 
at Zi is then given by 



The first factor in the numerator stems from the image vortex at the plane Y = with 
complex potential w(Z) = — z(K/27r) ln[Z — Zi] (which has to be present even without 
the circle), the second one is obtained by the circle theorem as the image of the 
original vortex at the circle. Finally, the potential of the remaining +K-circulation 




(27) 
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Fig. 1 The most simple nontrivial boundary problem solvable by the image technique: A 
unit circulation point vortex in the half space Y > 0, which is filled with liquid, moving near 
a (half-)circle. The boundary conditions are satisfied by an image vortex at the plane and 
two image vortices of opposite strength inside the circle. 



vortex inside the circle, contributing in the denominator of the logarithm, completes 
the image vortex system, again by the circle theorem. The first term in the de- 
nominator of the logarithm is incorporated into the static self energy of the vortex, 
Eseif = {mpoK? / in) In (\Zi — ^i|/C) , which is cut off by ^ and equal to half the energy 
of a vortex pair separated by \Zi — Zi\. The expression for the potential of (0) is thus 



iPc = -(«/27r) In (I (dV^i - Zi)/{dyZi - Zi)\) 



(28) 



which is just the counterpart of ipc in ( |2q ) in two dimensions. The complex counterpart 
of tpu, in turn, is 



Wu{Z) = u[ Z + — 



(29) 



giving the complex potential of flow from right to left (in the negative X-direction). 

The full energy of the point vortex, neglecting compressibility effects, and using 
the normalisation (21), is thence given by the expression 



E{Zi) ^ In 



(Zi - Zi) {dyZi - Zi) 



i{(P/Zi-Zi) 



Attu^ ( (f 

K \ Zi 



(30) 
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2.6.2 Conformal transformation 

A conformal transformation is a co-ordinate transformation leaving angles invariant, 
that is, the metric is multiplied by a conformal factor which is the transformation's 
Jacobian determinant. What we want to do is to map by a conformal transformation 
the region outside a boundary surface with varying curvature radius, lying in the 
(target) 2;-plane, to the domain outside the circle, which is in the (original) Z-plane. 
Any such transformation can be written as a holomorphic function of Z (save for 
singular points, such as vortex centers), in the form 

oo 

z^aoZ + Y.bnZ-'', (31) 

n=0 

where ao,6„ are some coefficients and Z = (iexp(ix) is on the circle (we omitted a 
constant, indicating a change of z-plane origin). 



2.6.3 Tie ellipse solution 

We would like to invert relation (^Tj) to obtain the solution for the boundary surface 
directly from that for the circle, which we have already obtained above. Easiest 
to perform is this inversion if we let 6o = 0, = for n > 1. Furthermore, if 
ao = 1 by proper normalizing choice of scale, we are led to the celebrated Joukowski 
transformation 



z = Z-f/4:Z, 



(32) 



which maps the outside of an ellipse with half-axes a, b (where a < 6) to the outside 
of a circle of radius d = {a + b)/2. The parameter I is defined by P = — a^. The 
inversion of this transformation is cast into the form 



Z = 



i (z + ^/z^ + P^ = ^ (sinhC + coshC) = ^ exp[C] 



(33) 



with the aid of elliptic co-ordinates, defined by (C = X + *^) 

z = I sinh C — I (sinh x cos rj + i cosh x sin 77) — x + iy . (34) 
The lines of constant 77 and x are confocal ellipses and hyperbolae, as follows from 



y 



P sin r] P cos^ rj 



1 



y 



P sinh x cos2 x 



1 



(35) 



The co-ordinate basis = d/dx, = d/drj is an ortho-basis with the conformal 
metric gtj = Z^(cosh^ x cos^ 77 -I- sinh^ x sin^ v)^ij ■ 

The normalized energy (|30| ) as a function of the elliptic vortex co-ordinates xi i Vi 
takes the form 



E{xi,Vi) = In 



b exp(xi - Xo)| sin 771 1 sinh(xi 



Xo) 



(sinh^(Xi - Xo) + sin^ 771)1/2 



Anuia + b) 
smh(xi 



Xoj|sm77i| , 



(36) 
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where xo = artanh(a/6) is the co-ordinate specifying the (half-)enipse surface at the 
boundary, i.e. the half-axes of the eUipse are given by a = ^sinhxoi b = ^cosh^o- 



For direct comparison with the half-ring case treated in subsection 2.5, which had 
a quite simple geometrical meaning, we will investigate the paths of constant energy 
E{xi,r]i) = Eq mainly in the case of small velocity u. The notion of 'small' requires 
some more care in two dimensions as compared to the cylindrically symmetric half- 
ellipsoid case in three dimensions. Whereas in the latter the velocity enhancement at 
the top is always 2u, in the 2d case it is 2{b/a)u. Hence it is required that we restrict 
the velocity at infinity to 

" « ^ , (37) 
such that the velocity at the top is well below the critical velocity 

^-^R' (38) 

which gives a measure of the onset of many-body quantum physics in the atomic 
superfluid helium II. For Fermi superfluids, the corresponding critical velocity is the 
pair-breaking velocity of Cooper pairs. We will henceforth refer to vl as 'Landau 
velocity' and generally scale velocities with it. Numerically, the actual Landau critical 
velocity of roton creation ~59m/satp~l bar equals vl if ^ ~ 2.7 A is taken.^ 

The approach of the vortex to the ellipse surface is closest at the top {rj = tt/2). Let 
us calculate the distance at the top as a function of the constant energy Eq. In general, 
a small distance interval ds is given by ds — \/ gxxixo) Sx = (a^ sin^ rj + P cos^ riy^^Sx- 
We define a quantity s by 

s = aSxc = aSxiil = 7r/2) • (39) 

At the ellipse top, s — Ss{r/ — 7r/2), provided that Sxc = Xiiv = ^/2) ^ Xo ^ Xo- 
The actual distance of the vortex to the top is somewhat different if this inequality 
does not hold. We do not further dwell on this difference here for the sake of simple 
argument and take (^) as a definition of the quantity of distance s. 

Evaluating ( p6| ) at 77 = 7r/2 for small velocities, and assuming that ^0 is small 
enough for Sxc ^ 1 to hold, we get Sxc — Cexp[£'o]/(a + b) and therefore 

s ^ exp[£;o] S.-^ exp[Eo] ^ (77 = 7r/2, u smah) . (40) 

On the other hand, far away from the half-ellipse, at the flat boundary, the distance 
is, in the low velocity limit, given by Ss ~ (1/2)^ exp[_Bo]. We observe that the vortex 
comes closer to the ellipse by a ratio 2a /b (for 6 ^ a) than it was far away from the 
ellipse. If a 6, that is, in the case of the circle, the distance remains the same. 



'^It should, however, be pointed out that the very notion of (classical) velocity becomes questionable 
on ^-scales. On these scales, it is more appropriate to refer to (mean) current densities. 
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2.6.4 Geometric restrictions 

Now, if we are bound to remain in the realm of the hydrodynamic description of a 
collective co-ordinate vortex, we have to impose that the total distance As of the 
vortex to the ellipse always exceeds a quantity O(^): 

As= r'6s>O{0- (41) 

Xa 

The quantity O(^) means 'a value in the order of ^ by definition', as there can obvi- 
ously be no sharp distinction between the 'inside' and 'outside' of a quantum vortex. 
Additionally, the choice for the lower limit value of the total distance / 6s in units 
of ^ depends on the value of the core constant Co — 0{1), parameterizing the many- 
body core structure in the vortex energy logarithm ln[i?/(^exp[Co])]. For the point 
vortex considered, which has R ~ 2y,expCo = 1; for a ring vortex we chose the 
parameterization 8exp[Co] = exp[C], cf. p^ . 

It is apparent from (|40| ) that if the energy of the vortex is sufficiently small (in 
particular, if it is zero) and the ratio b/a is sufficiently large, the condition (^ ) will be 
violated and a hydrodynamic collective co-ordinate description invalid. Thus, there is 
a minimum vortex energy required for the whole formalism we employ here to retain 
its validity. Expressing this energy as a function of the parameters involved, we have, 
for general velocities {Sxc = Sxiv — ''■/2)): 



Eo = In 



— expSxc tanh(5xc 



smh 5xc (42) 



K 



~ In 











2u f b\ s 

— 1 + - - . 43 

the last line valid provided that we can approximate, with reasonable accuracy, e.g, 
tanh^Xc — 5xc- The energy Eq is the energy needed by the vortex to remain com- 
pletely describable in terms of a collective co-ordinate on its way along the ellipse. 
The quantity s is defined in (^9|), and has a lower bound related to the prescription 
(©■ 

We have depicted the normalised potentials corresponding to different velocities in 
units of and ratios s/ a in Figs. and ^j. The potentials in these figures correspond 
to the real space ellipse shown in |5| which has b/a ~ 5.7 and where s/a « 2. For 
clarity, we have additionally displayed the shape of the barriers in the direction of the 
y-axis (i.e. the e^-axis at 77 = 7r/2) in Fig. IJ. Whereas for u — 0.04 the semiclassical 
approximation is still applicable, for u = 0.08 this is no longer the case. It is approx- 
imately for this velocity that the barrier tends to zero and the tunnelling distance 
approaches ^. One can see here explicitly that indeed represents a sensible measure 
for a critical velocity, because the local velocity at the top is ~ ul, if the velocity at 
infinity takes the value u ~ 0.088. 



The situation we encounter in real space is shown in Figure |5|. We visualize the 
quantum core size as the shaded area around the vortex center, having the 0(0" 
radius used in (^l]). The vortex on path 1 is not able to pass the ellipse without some 
part of this shaded area covering the ellipse, but the vortex on path 2 with energy Eq 



18 



Ann. Phys. (Leipzig) 9 (2000) 1 




Fig. 2 Shape of the potential barrier ( |3q ) with the choice xo = 0.175 for the velocity 
u = 0.04 in units of vl ~ k/2ii£^. The corresponding real space ellipse with h/a ~ 5.7 is 
shown in Fig. ^. The ratios a/^ = 2, s/^ = 1 and thus s/a = 1/2. The zero of this 
normalised potential energy is shifted by Eo, defined in (^2|). The contour lines of the vortex 
paths of constant energy shown are those for E — Eq. 



avoids the ellipse surface completely. That the vahdity of the hydrodynamic approach 
enforces that we introduce another geometrical quantity, s, the vortex distance of 
closest approach, which implies a lower bound of vortex energy, Eo{s), is an observation 
of general character. It is of relevance for any attempt to describe tunnelling in a 
realistic, non-spherical geometry, i.e. when the boundary and thus the path of the 
tunnelling object near it is not of 5" symmetry. A hydrodynamic collective co-ordinate 
description is valid only if the quantum core structure of the tunnelling object is not 
touched upon during its motion along the boundary. A pinning potential for the 
vortex moving in the superfluid stems in general from some flow obstacle, in our case 
the ellipse. Any phenomenological ansatz for a pinning potential usually employed in 
tunnelling calculations, which has curvature perpendicular to the applied flow larger 
than parallel to the flow will have to take into account that the object can approach 
the surface within its core size. 

We now come back to the determination of the two paths of constant energy we 
need to construct the closed path in Euclidean phase space. Assuming (5^ 3> 1, we 
get the path far away from the ellipse to which the vortex has to tunnel. It is, with 
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Fig. 3 Same potential as in Figure EL but for velocity u — 0.08. 



decreasing magnitude of u, given by the relations (as a reminder, 6 3> a): 



^exp[xi]|sin77i| 



■In 



27ru 

^In 

u 

^In 



2Yi 



■exp[xi]|sinryi| 



VL 



u exp[£'o] 



(44) 
(45) 
(46) 



The relevant path is the solution which has Yi 6 (the second solution of the equations 
above is just the vortex trapped at the boundary, far from the ellipse). The last relation 
is exactly analogous to (p3|), which is valid in the low velocity limit of three dimensions. 

The path at the ellipse with energy Eq does have approximately constant (5x <C 1 
if s <C a. The shape of this path can be obtained (to lowest order in 5x) by solving 
the equation 



(a + 6) I sin ?7i | = 



In 



(a + b)5x 
^exp[i?o] 



(47) 



Around the ellipse top, within a large range of ry-values, the solution is approximately 
given by 5xc — s/a. This can also be inferred from Figure |[ 
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Fig. 4 The potential barriers of Figs. |2| and ^ in the predominant direction of tunneUing 
along the e^-direction at 77 = 7r/2 (the j/-axis). Left: u — 0.04, right: u — 0.08. 



2.6.5 The tunnelling area 

To describe the motion of the vortex in momentum space, we choose a gauge for the 
momentum which is most appropriate to the symmetry of our problem. In Carte- 
sian co-ordinates, this is Px — hpgY, i.e. the gauge momentum equals the physical 
momentum (i.e. Kelvin momentum) for an unconstrained vortex in the bulk super- 
fluid. We will evaluate the action in the limit of small velocities. In this limit the 
path to which the vortex has to tunnel (named in what follows Iat), is given by a 
constant, Yq in (|4^). We have also seen at the end of the preceding section that 
the vortex remains approximately on the same ellipse, having the elliptic co-ordinate 
Xi — Xo + <^Xc(7''/2) — Xo + s/a, while it is moving around the ellipse (path 2 in Figure 
H). Calling this path Ye, we can write 

Yi = Pcosh^{xQ + Sx) + tanh-^ixo + Sx)K^ 

~ P cosh:^ {Xo + Sxc) + t&nh-^ {Xo + Sxc) , (48) 
using the imaginary co-ordinate K — —iX. Then, the integral in (|^) takes the form 

oinc 

-f- = 2po iY^-YE)dK, (49) 

where Km is the point at which the trajectories Y}v and Ye meet in complex phase 
space. This point, determined in the low velocity limit by the solution of Ye = Yq, 
can be shown to be 



Km ~ tanh(xo + 5xc) Yq 
'a ^ s 
- b a 



)^ln 




/ u 


.M exp[£'o] . 



(50) 



neglecting the first term in (p8|), which is possible if Yq '3> ?cosh(xo + Sx) = &x- 
last line is valid provided that both xo, Sxc ^ 1- 



Now, the integral (49), with = a(xo + Sx) = 'sinh(xo + Sx), takes the form 



cine / u 

-j- ~ 2po [KmYo -^Kl]^po tanh(xo + Sxc) Y^^ . (51) 
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Fig. 5 Two vortex paths of constant energy near the ellipse. Whereas the vortex on path 
1 with approximately zero energy, iJ ~ 0, does not manage to pass by without coming closer 
than O(^), the second one, having energy E = Eo, defined in (p^, is able to do so. The 
velocity u is to be understood that of the flow 'at infinity'. 

The above integral corresponds to the shaded area in Figure ^. Neglecting the small 
cusp at the bottom reducing slightly this area, as we did in the third line of (|l|) , leads 
to a volume lav^ of the form 

^ p,n^^^ ^ P,PV^^') . (52) 

Here, the effective sharpness /? = tanh(xo + 5xc) is a measure of the maximal ec- 
centricity e ^Jl — fi"^ ~ 1 — (l/2)/9^ the vortex path is allowed to have, under the 
condition that the vortex should remain completely within the hydrodynamic, collec- 
tive co-ordinate domain. 

The area of Fig. || thus does have a lower limit. In general, (3 characterizes the 
effective dimension of the vortex escape path, i.e. the relative degree to which this 
path is confined to n dimensions by the presence of an asperity which is effectively 
71-dimensional. The tunneling volume is that for a vortex escape path of 0(rf — 1) 

®We refer to the contribution of S]P-'^ as a volume contribution in general, according to (|l7|). In 
two dimensions, of course, a volume is an area and an area is a length in conventional terms. 



22 



Ann. Phys. (Leipzig) 9 (2000) 1 




K=-iX 



K = \ 



Fig. 6 The closed vortex path giving the action (^l|) in the low velocity limit. The first 
part Ye corresponds to the analytically continued path 2 of Fig. |^ along the ellipse surface. 
The second part of the closed path, Ym — Yq, represents the border line to a free vortex. For 
higher velocities, this simple triangle shape is deformed. 



symmetry, which is the highest possible symmetry if one preferred direction, namely 
that of the external current, is given. In the d = 2 case treated here at length, we have 
n >1, vjf'' — Yq for a single vortex and correspondingly V^^'' = 2Yq for a vortex pair. 



In three dimensions {d = 3) we have V^^' = (27r/3)i?g for a half-ring with radius Rq 
and double this value for a full ring. The effective sharpness will be reduced {i.e. the 
value of P larger) for the 3d case (cf. the equation (25)), in analogy to our analytical 
findings for the vortex in the plane. In order of magnitude, the sharpness will in any 
dimension be given by the product of the ratios of the curvature radii parallel and 
perpendicular to the flow of the allowed vortex path. 

We conclude with an estimation of the kinetic contribution in the action, using 
the last line in (^5|). We will also use that in lowest order of perturbation theory the 
velocity of the vortex equals the local flow velocity Vg in the incompressible superfluid 
(the Kelvin-Helmholtz theorem). The velocity maximum (at the ellipse top) we set 
t^max = max(lwJ). Then, 



5f"(i?) - hpo^^^^i{Qa/Cs)dX'' 



< 2hpQ 



2tt 
ln(-. 



2tt 



J. ^inax 

■4 Kj, 
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= Mlli (53) 

TT Yo Cg 

The contribution of the vortex mass is thus at least suppressed by two small factors: 
^/Yq and Vmax/cs, contributing both in equal measure to the fact that S'g'" is negligible 
as compared to S^^^"^. 



2.6.6 Summary 

We now summarize this extended investigation of the analytically soluble 2d problem, 
emphasising in particular its crucial outcome. 

Starting from the 'electrostatic' problem of a single point vortex situated near a 
half-circle at a boundary, we derived the vortex energy by conformally transforming 
via the (inverse) Joukowski transformation to the half-ellipse solution. We concluded 
from the general expression for the vortex energy (|3^) , that it is necessary to introduce 
the geometric constraint ( p] ) , expressing the limits of the hydrodynamic collective co- 
ordinate formalism under consideration. Calculating the tunnelling volume (area), we 
have seen that it assumes the general form (52). The tunnelling volume (area) fi'-^-' 
cannot be reduced to a lower bound given by the sharpness of the ellipse, tanh^Oi but 
has to be larger, P = tanh(xo + ^Xc), if we are bound to remain within the domain of 
the approach which has been employed. 

It is to be noted that we expressed the energy of the vortex in units of the energy 
h'^ Pq / {4:Trm) = mpoK'^/iir ~ 0.82 K/A (at p ~ 1 bar and in three dimensions). For 
realistic values of the parameters in (^^, the values of Eq cover the same range as 
the phonon-maxon-roton spectrum. Energywise, the trapped small scale vortex thus 
cannot be distinguished from an elementary excitation of the superfluid. It could have 
been excited thermally and remained trapped at a pinning center during the cool-down 
of the superfluid to very low temperatures. 

From the above analysis it thus follows that it is semantically in general not quite 
appropriate to employ the widely used term 'vortex nucleation' for vortex tunnelling 
investigated within a hydrodynamic, collective co-ordinate theory. If we define 'nucle- 
ation' to mean creation from the zero of energy, we have seen that nucleation is not 
amenable to such a description under general circumstances. Experimentally, it will be 
impossible to distinguish the tunnelling of a small energy vortex at a rough boundary 
from the true nucleation event of a nascent vortex there, if no direct means to control 
the microscopic dynamics can be provided. It remains, of course, to be explained how 
a vortex of such small energy and size (i.e. one with small distance of O(^) to the 
boundary), can be defined and described quantitatively. The present approach can 
only fix the tunnelling rates of a vortex if one presupposes that the vortex can be de- 
scribed by a collective co-ordinate along its whole path. This, however, will be true if 
the size of the vortex is sufficiently larger than O(^). Then, the dissipation-free motion 
of the vortex with constant energy and the purely geometric nature of the problem 
lead to strict bounds for the possible values of the semiclassical tunnelling exponent, 
directly related to geometrical quantities. 
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2. 7 The prefactor 

Any discussion of quantum tunnelling is incomplete without at least an estimation 
of the prefactor A{E) in (|l]). If the quantum fluctuations, or better, the indetermi- 
nacies of the vortex position and momentum vanish, the tunnelling probability does 
the same, because the very process of quantum tunnelling stems from this quantum 
uncertainty of the vortex in phase space. The crucial advantage of our WKB-like in- 
vestigation lies in the fact that the behaviour of the exponent, the Euclidean action of 
the instanton, dominates any dependencies of the prefactor on observable quantities 
as long as Se{E)/h ^ 1, i.e. as long as we are in the semiclassical limit. We will write 
(^ in the form 

P{E) =exp[-SeiE)/h + liiA{E)] , (54) 

to compare ~Se{E)/h with the dimensionless quantity \nA{E) (the prefactor A{E) is 
originally in units of Hz). 

2. 7. 1 Estimations 

Apart from the considerable difficulties in evaluating prefactors in general, an accurate 
calculation of A in a dense supcrfluid like helium II is in principle not possible at 
present, due to the lack of a microscopic theory. It is, however, feasible to get an idea 
about the value of this prefactor within about two orders of magnitude. 

Let us begin with a physical picture for the origin of the prefactor in the semiclas- 
sical limit. The simplest possible idea about the prefactor is gained by considering 
the frequency Ua of a particle oscillating in a metastable well. Then, within about 
one order of magnitude, A ~ In the thermal activation limit, i.e., in the 

Arrhcnius law case we have more exactly P = {uja/'i.'K) exp[— [/ /fc^T^], where LOa is the 
frequency of oscillations at the metastable well bottom against a barrier of height U . 
The frequency Va = ^a/2.'K can generally be understood as a measure of the number of 
times per second the vortex bounces against the potential barrier, trying to get free. 

We have no possibility to describe the vortex state (at the boundary) quantitatively, 
but we are able to conclude on the order of magnitude of uoa, if we take into account 
that there exists a surface layer of vorticity, of width ^: Because the superfluid density 
goes to zero at the boundary and heals back within ^, the energy needed for the 
activation of vortices vanishes within this distance |2^, The frequency of motion 
of these vortices should then be of order 

Wo = 4.87 • 10"sec-ir^[c^Lf] (helium II), (55) 

which is the cyclotron frequency of vortex motion [3C| ] . We scaled ^ with the Lennard- 
Jones parameter ctlj = 2.556 A of the ^e atomic interaction. The frequency loq is 
the natural vortex frequency associated with the scale ^ alone. 

If we wish to connect Wj, with the quantity relating the strength of interatomic 
forces and compressibility, namely the speed of sound Cg , which is involved in the vortex 
kinetic term, the following phenomenological treatment is of use. The expression for 
the vortex self energy is logarithmically divergent. If we consider ^-scales, this may 
be cured in a heuristic manner by regularizing the logarithm under the condition that 
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the vortex energy be zero at the boundary (that is, at F = 0) This leads to the 
following total energy of a point vortex in the laboratory frame, using the normalisation 

E{Y, X) = 1 [(X/csf + In (1 + (y/0') - SttuF/k) . (56) 
The minima of the potential, existing for low enough velocity, are situated at 

l^n^in = ^ (l - VI-^Hu/vl)') - f-V = 2e - , (57) 

where in the final result we used an approximation for u <^ vl- The curvature of 
the potential is ^~^(1 — (^/0^)/(l + (^/O^) ^^d hence at the minimum, for u <C wl, 
approximately 1/^^. We have thus found that the spring constant of an 'elastic object' 
vortex should scale with a quantity of order 1/^^. Because the mass is in the 
same units, the frequency of oscillation is therefore 

c..^^^9.4.10"Hz^i|i^ (58) 

The ratio of the two estimates (p5|), ( p8| ) is given by wq/ws = {n/ Cs) / {iri). They thus 
coincide in order of magnitude because the relation (p^, valid in helium II, holds. It is 
quite obvious that both estimates can only give a rough idea about the actual attempt 
frequency in the dense superfluid, whereas at least the estimate (|5q) s hould be quite 
accurate in dilute superfluids (for 'nonrelativistic' vortex velocities [|3l|]). However, by 
argument of continuity, we do not expect the analytical relation for the frequency oja 
as a function of the parameters ^, Cs to make abrupt changes, if we increase the density 
of the superfluid. We suspect that the density increase does not cause a change of more 
than, say, one order of magnitude in uja, than that predicted by the above estimates. 

Now, relying on the estimate (|5^ ) and scaling the frequency with 1.13 • 10^^ Hz (the 
roton frequency at p ~ 1 bar), so that P[E) = F(_E)[Hz] and [cs/C}n = (cs/C) [1-13 • 
10"'^^ Hz], by use of ([l^), we arrive at a tunnelling probability having the appearance 
of (H): 

P[E) = exp {-po(2^f^^'')+ln(---)CS('')) + 27.75 + ln[(cs/On]} 

~ exp |-po PY;^ + 21n(- • •) |- ^\ + 27.75 + ln[(cs/0„] 



exp 1-2^7^^) ( 1 + ^ + 27.75 + 

I V TT Yq Cs J 



ln[(c,/0«] , 



(59) 



where in the second line we inserted the results of the 2d half-ellipse problem in 
equations (^2|) and (p3|), and the third line employs the Bohr-Sommerfeld quantization 
of (|19|). 
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3 Comparison to experimental results 

Having determined wliat we can expect for the magnitudes of tunnelling probabilities 
we now come to discuss how experimental arrangements to measure such tunnelling 
events are constructed, and what the findings of such experiments are. In most of the 
measurements (cf. Q~[Q), the apparatus to measure phase-slip events is structured 
like that shown in Figure ^ We note in passing that such an apparatus with an 
effective torus geometry can be used to measure external rotation {e.g., that of the 
Earth) quite sensitively (see Refs. pl[-[5^). 



md,r]d,k,F 



Sk, h 



He II 




Fig . 7 The schematic arrangement of the ZAV (Zimmermann-Avenel-Varoquaux) oscillator 
|44 |. A membrane, generating a current It by an (electrical) driving force F, has mechanical 
parameters mass rud, stiffness k and damping constant rjd. This current divides into a current 
Iw = It{Lk/{Lw+Lk)) through a micrometer-sized small hole and a current Ik ~ It{Lw/{Lw + 
I/fe)) through a long channel. Here, Lu),fc = lw,k/{'rnpoSw^k) are the hydrodynamic kinetic 
inductances, with l^^k the lengths and Sm,k the cross-section areas of the micro-orifice and the 
long channel, respectively. The circulation threading the hole and channel, Iq, is quantized. 
The whole cavity is filled with He II. 

Essentially, it is observed in these experiments that at a well-defined value of the 
amplitude of the diaphragm, driving the current through the micro-orifice (the small 
hole on the left-hand side of Figure 0) , there is an instantaneous (on the scale of the 
driving frequency) breakdown of the diaphragm amplitude, which is quantized. This 
quantized dissipation event is associated with a vortex generated at the orifice walls, 
subsequently crossing all the streamlines of the fiow, thereby causing a phase slip event, 
which draws a quantized amount of energy from the fiow. We will first discuss these 
phase slips jSSt , which give the crucial physical argument for the interpretation of the 
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experimental results. 



3.1 Phase slips 

The picture of phase slips is based on the fact that the phase and particle number are 
canonically conjugate^ that is, 

[N, 0]=i. (60) 

The quantum mechanical equation of motion of the phase is hence given by 

Taking the thermodynamic average of this equation, we see that the time rate of 
change of the phase equals the (negative) local chemical potential (defined hy p, = 
fi + {l/2)mv'^, where ^ is the chemical potential in the superfluid rest frame), divided 
by h, a relation which became popular under the name 'Josephson- Anderson equation' 
[ p4| , |5^ . We now consider the time average of this thermodynamic average over a long 
time span t and take the difference of the results for two points A,B in the superfluid: 

-(A(r^)-M(r-B)>, = (hf'^^^^ 



\ dt 



= lim 

r — i-oo 



!^ fdt 4 ( / ve-ds 



T 



'^^ \"'Cab 



The last relation tells us that the the negative chemical potential difference between 
the points A and B, divided by Planck's quantum of action, is equal to the number of 
vortices crossing the line Cab, joining A and B, per unit time, dn^/dt. A single phase 
slip process, caused by one migrating vortex, may be visualized as shown in Figs. 
^ and ^. The Figs. [H [o] represent a pronounced simplification of the actual vortex 
motion process. We visualize in these pictures a point vortex moving on a straight line 
across the superflow through the orifice, which gives a highly symmetric view of the 
process. The real process of vortex half-ring motion is investigated in Refs. j4^, 
The principal global (topological) features of importance however remain untouched 
and independent of the actual vortex motion trajectory: The phase difference between 
two stationary states (times i <C ii and t ^ in Figs. ^ and ^ is exactly 27r and 
the process always sucks the same amount of quantized energy from the flow, given 
by Ai? = mponSwVc — nJc |^6|| , with Vc being the (mean) critical velocity of flow 
through the micro-orifice, at which the vortex migration process sets in, and Jc the 
corresponding mass current (see for a derivation below). 

In Figures || and ^ we represent a quantized vortex crossing the micro-orifice 
designated with quantities S^, Iw in Figure ^ The lines emanating from the vortex 

'This is true in the hydrodynamic limit, after averaging over a cell much larger than the atomic 
size. In the microscopic domain, we encounter consistency problems related to the general problem 
of the existence of quantum mechanical phase operators (see, in particular, js?! ; also js^). 
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Fig. 8 The early and intermediate stage of a phase shp process in the orifice |55| 



center (black dot) in Figures ^ and ^ are lines of constant phase, standing perpendicular 
on the orifice wall (which is a streamline). For ease of representation we have chosen 
the branch cut of the phase to be exactly parallel to the direction of motion of the 
vortex. The whole phase slip process is then maximally symmetric. The shaded areas 
represent the walls of the orifice Sw , Iw on the left hand side of Figure 0. The points 
A, B are chosen to lie sufficiently far away from the orifice, as indicated by the dots. 
Initially, at t = ti (first drawing), when the vortex starts on the left side of the orifice, 
the phase difference between A and B, Oj^ — Bg, is zero. When the vortex reaches the 
line joining A and B, the phase difference is tt (second picture in Figure Finally, as 
shown above, the vortex has crossed all the streamlines through the orifice, flowing, as 
indicated in the first Figure, in the vertical direction, and disappears to the right. The 
drawings above give a representation of the effect of migration of the order parameter 
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Fig. 9 Final stage of the phase slip process. 



zero, the vortex, across the matter flow through the orifice, by means of the phase 
of the order parameter. The physical result of this migration is invariantly given by 
the energy change, expressed in terms of the critical mass current for triggering 
the generation of the vortex, Ai5 = kJc, and does not depend on a representation 
in terms of the phase (or, for that matter, on the choice of the branch cut), i.e. is 
gauge invariant; it is also invariant under changes of the location of the points A 
and B (as long as they are situated far away from the micro-orifice). Multiplying 
the thermodynamic Josephson relation with the critical number current through the 
orifice, results in a time rate of energy decrease of the external flow driven by the 
oscillating membrane in Figure |^, -E = (KJc/27r)^, from which the energy change 
for one phase slip, necessarily causing 6 to change by exactly 27r, follows by time 
integration.^ 



3.2 Principal findings 

In the experiments, the critical velocity Vc — Vc{T) as a function of temperature T 
is measured. The critical amplitude of the diaphragm, corresponding to Wc, is that 
for which there is a diminished amplitude in the next half-cycle of oscillation. The 
Figure |l^ shows a typical experimental run of measured resonator amplitudes. An 
important feature of Vc is that it has a statistical distribution, which has also been 
recorded. The results, obtained from the statistical analysis of the series of phase slips 

^The resulting expression for the energy change neglects a small correction due to the diminishment 
of the external current during the phase slip feq. 
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Fig. 10 Elementary phase slip processes in two different runs The resonator amplitude 
of the diaphragm, which is ramped up in time by applying a voltage bias on the coated 
diaphragm, is recorded every half-cycle (triangles pointing up and down, respectively). A 
phase slip occurs if there is a breakdown of the amplitude to the next half-cycle. 



like those in Fig. are shown in Figs.|Tj, 12, ^ (received from E. Varoquaux and 
reproduced here with kind permission). The sahent results are that the mean critical 
velocity first rises linearly with temperature and then saturates at Tq ~ 150 mK. Cor- 
respondingly, the statistical width decreases linearly and saturates at approximately 
the same temperature. 



3.3 Interpretation 

According to the phase slip picture we have developed above, a possible interpre- 
tation of the experimental data is as follows. A vortex half-ring, standing with its 



axis antiparallel to the flow |28|, is generated at the orifice wall and expands under 
the influence of the diverging flow fleld through the orifice. During this process, it 
crosses all the streamlines through the orifice, completing the phase slip, and is finally 
transported away from the orifice by the flow. The fact that there is a certain criti- 
cal amplitude of the diaphragm for this procedure to happen, can be associated with 
the fact that there is a potential barrier opposing the process. Furthermore, the fact 
that the critical amplitude {viz. the critical velocity Vc), has a statistical distribution, 
supports the idea that the existence of a phase slip critical velocity has the statistical 
origin of barrier crossing events. Additional support is provided by the linearity of the 
amplitude respectively its distribution with temperature, a signature of thermal acti- 
vation over barriers ||l^ . It is also measured that the (average) flow velocity through 
the orifice as a function of temperature, needed to trigger the phase slips, as well 
as its statistical distribution, saturates at a temperature of ~ 150 mK. This can be 
ascribed to a non-thermal process of surmounting the existing barrier: The possible 
explanation of the observed behaviour is quantum tunnelling of half-ring vortices at 
boundaries. 

Within a phenomenological approach , a model of half-rings with axis an- 

tiparallel to the applied flow, and standing perpendicular to the walls of the orifice, 
has been developed. The Hamiltonian is essentially that in (|5^), save for the kinetic 
term. It turns out that, to make this model conform with the available measurements, 
it is necessary to postulate, a) that the coherence length increases to ^ ~ 9 • • • 10 A at 
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Fig. 11 The critical velocity through the orifice (in terms of instrument units of the critical 
amplitude of the diaphragm) |^^. The line is a fit to a (thermal) 'half-ring model', discussed 
in the text. The temperature region below 200 mK is magnified in Fig. 



boundaries and b) that the vortex half-ring be described by a collective co-ordinate on 
0(^)-scales. This approach must consequently be understood as the parameter-fitting 
of a simplified model to the available data. From a more fundamental point of view, 
it does not describe crucial, indeed salient features of the actual problem: 

i. The model does not incorporate Galilean invariance violation. We have seen 
that this is a necessary requisite for any (hydrodynamic) formalism making, in 
particular, use of the notion of velocity, to describe the quantum tunnelling of 
vortices at temperatures close to absolute zero. 

a. It cannot be reasonably expected from a description of the entity vortex on ^- 
scales to make sense for the tunnelling exponent beyond crude order of magnitude 
estimates, like in the case of the prefactor (where these estimates are sufRcient). 
For the semiclassical tunnelling exponent, all nontrivial dependence on coherence 
length (respectively microscopic) physics, in whichever form, should be excluded, 
such a dependence only being allowed in the form of an ultraviolet cutoff. 

in. It is certainly not permissible to neglect the dynamic influence of the kinetic 
energy of the vortex in the tunnelling exponent, if we approach scales of order 
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Fig. 12 The critical velocity statistical distribution. 



the coherence length. 

Let us now further analyse the experimental outcome. First of all, we rely on the 
hydrodynamic relation (p3| ) to deduce the scale of materialisation of the nascent vortex 
half-ring. Scaling the velocity with 10 m/sec, the order of magnitude of the (local) 
flow velocity, we get 

1.59 nm / 9.89 \ 

" ~ M[10m/sec] U[f^Lj]u[10m/sec] j ' ^ ' 

where use was made of the Roberts-Grant result [M for C = 1.615, valid within Gross- 



Pitaevskii theory [^-|63 . In the experiments, the measured velocity through the hole 
is 5-10 m/sec. This can only be measured as an average over the cross-section of the 
orifice, locally the velocity can of course be higher. Nevertheless, it can be concluded 
that the radius Rq should be of the order of nanometers. This small mesoscopic scale 
makes it difficult to decide if a hydrodynamic formalism is applicable in a rigorous 
sense (one should also bear in mind that the formula above is strictly valid only in 
the low- velocity limit). In particular, the value of the coherence length is not exactly 
known under the circumstances considered. The neglect of any large density variations 
{Sp/po — 0(1)) in the formalism makes it necessary at least to assume that Rq ^ ^, 
so that the knowledge of ^ is crucial indeed. 



Next we consider the value of the cross-over temperature Tq ~ 150 mK. (ReL 
19[ reports a value Tq ~ 200 mK.) The crossover temperature between thermal and 
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Fig. 13 Magnified portion of the temperature region below ~ 200 mK [^s],^]. The figure 
shows the infiuence of minute impurity concentrations of ^He, given in ppb, on the critical 
velocity of vortex phase slip for very low temperatures (for clarity, the curves are shifted 
downward with respect to one another by 200 instrument units). The dashed curves are a 
fit to the data, according to a phenomenological theory that "'ffe atoms, bound to the vortex 
core at the tunnelling site, reduce the critical velocity [^. The dash-dotted curves are the 
high-temperature linear fits shown in Fig. For the purest sample (0.9 ppb ''He), the 

critical velocity is temperature independent within experimental accuracy down to the lowest 
temperature ~ 15 mK. 



quantum behaviour gives in general a measure of the equality of quantum-mechanical 
and thermal energies of the 'particle', trying to surmount the barrier with the aid of 
these energies [ [73| . For zero damping (dissipation) the crossover temperature is given 
by hujh = 2'!TkBTQ, where LUf, is the frequency of oscillation at the top of the barrier, 
connected to the trivial solution of the Euclidean equations of motion of the 'particle' 
sitting at the top (the bottom of the inverted potential) This leads to 

W6 = 8.2-10^"Hz-To[100mK]. (64) 

If we compare the value of ui, in ( |6^ ) with the frequencies ( p5| ) , (^8|) , we see that it is 
smaller by a relatively large factor, up to one order of magnitude, provided we assume 
the coherence length to have its bulk value, which is about ctlj- However, a direct 
comparison of the experimental value and these estimates cannot give more than an 
order of magnitude agreement. This has several reasons. First of all, these estimates 
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can give only a very approximate idea about the true dynamical behaviour of a quan- 
tum many-body vortex near the boundary. It is conceivable, for example, that the 
effective 'spring constant' of the vortex against deformations is lowered compared to 
the semiclassical estimate in ( |55| ) , because of the many-body quantum uncertainty of 
its position. Then, the prefactor is in general a function of driving velocity u and 
temperature T p5[ |, and at the measured crossover temperature in terms of the criti- 
cal velocity not necessarily equal to its value at zero temperature. Finally, as already 
mentioned, we do not know the (effective) value of the coherence length at the bound- 
ary. It is conceivable that the value of ^ is enlarged as compared to the bulk, because 
of boundary conditions, i.e. the depleted superfluid density |Q. What one can thus 



definitely claim to have observed from Fig. 13, is that there exists a crossover tem- 
perature from a temperature dependent to a temperature independent regime, whose 
energy equivalent fc^To is in order of magnitude agreement with the estimates for the 
quantum oscillator energies hw^ and hws. To conclude this section, we give an idea 
about the number of particles involved in the tunnelling procedure, therein following 
the statement of equation (^9|) and the scaled tunnelling probability in eq. ( ^ ) . As- 
suming from the above discussion that the prefactor can vary in its order of magnitude 
between A w lO^*' • • • 10^^ Hz, its logarithm InA ~ 23 • • • 28. Corresponding to these 
conceivable values of the prefactor, the total number of particles in the tunnelling vol- 
ume should then be somewhere in the range 

Ar(3) = 4 - • - 6, say, for tunnelling events 
to be observable within a reasonable span of experimental time. Again, like the value 
for the materialisation radius Rq in ( |63| ) , this number indicates a rather small scale of 
tunnelling. 



3.4 Concluding evaluation of the experiment 

We have seen that the available data on critical velocities of phase slips can be inter- 
preted to be in phenomenological accordance with the picture of the quantum tun- 
nelling of vortices at boundaries below some temperature Tq. But the very fact of 
tunnelling at boundaries certainly needs further proof, so that the predictions of tun- 
nelling theory can be compared to that of classical instability mechanisms, which can 
be temperature independent as well. Theoretical investigations in this classical di- 
rection are found under Refs. |6^-|0 (vortex nucleation as a process of classical 
flow instability in experiments using 'Tfe-B was discussed in Ref. [|ll[). Provided that 
the hydrodynamic, large scale picture we developed here is applicable with sufficient 
predictive power for the actual materialisation scales of the vortex, one such proof 
could consist in the comparison of critical velocities for chemically identical orifices 
of equal global sizes, having different (microscopic respectively mesoscopic) surface 
structures. If the result of such measurements is negative, i.e. there is no reproducible 
difference in critical velocities, there is no quantum process taking place describable 
by hydrodynamic means of semiclassical tunnelling at irregular boundaries. 
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4 Aspects of vortex tunnelling in Fermi superfluids 

The analysis so far has been concerned with dense superfluids which are uncharged, 
and in which the fundamental constituents, i.e. the particles carrying the superfluid 
current, bear no internal degrees of freedom like spin. The spinless elementary bosons, 
which are '^He atoms, form the only hitherto known example of such a superfluid. 
There exists, however, a large number of superfluids which are constituted by elemen- 
tary fermionic particles. By far the most of these superfluids are charged: The charge 
carriers in superconductors represent a charged superfluid. Besides He II, the only 
other charge neutral dense superfluid known in laboratories on Earth is its isotope 
■^e. In the following, we give an overview of some general features of these Fermi 
superfluids, with particular emphasis on the vortex dynamical equations. The main 
intention of this section is to set the complications arising in the treatment of vor- 
tex motion in Fermi superfluids along the comparatively elementary hydrodynamic 
problem of unpaired bosons in He II. 



4. 1 Introduction 



For all Fermi su 
Cooper pairs |71 



Dcrfluids, there has to exist a mechanism binding the fermions into 
, which constitute the bosons of the superfluid condensate. In dis- 
tinction from the elementary *He bosons, these effective particles, arising from paired 
fermions, have in general the internal degrees of freedom spin and angular momentum. 
According to the symmetry of the order parameter, these superfluids arc classified to 
be S-, p- or d-wave superfluids. In the case of isotropic superfluids, the value of the 
internal angular momentum L = 0,1,2 of the Cooper pairs corresponds to s-, p- or 
d-wave, respectively. A relatively weak effective interaction binding the fermions to- 
gether leads to coherence lengths ^ 3> kp^, where kp signifies the Fermi momentum. 
The corresponding Cooper-pairs are then objects bound together over distances by 
far exceeding the microscopic scales relevant for many-body quantum mechanics. One 
consequence is that, e.g., the Ginzburg-Landau or mean-field levels of description, not 
suflicient to describe scales of order ^ in He II, are indeed useful for such paired Fermi 
superfluids on these scales |^]. The notable exception are high- Tc superconductors, 
where kp^ > 0(1). More properly, this is to be written as ^ > d^, where cIl is the 
lattice spacing, because the Fermi surface is in general not a single continous surface 
for the high- Tc materials. 

The fact that the Cooper pairs have internal angular momentum and/or spin leads 
to a richly structured order parameter, which can support symmetries much more 
complex than the global and local U(l) symmetries associated with particle conser- 
vation and the electromagnetic field. In particular, as one of the most important 
consequences of these enriched symmetries, there can occur nodes of the energy gap in 
momentum space for quasiparticle excitations above the superfluid ground state 173]. 



4.2 Vortex motion in Fermi superEuids 

Around a vortex line, there exists a potential well for quasiparticles, i.e. the pair 
potential is position dependent, A = A{x) (we designate Aoo to be the bulk value 
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of the gap, infinitely far away from the vortex hne). This leads to the existence of 
bound quasi-particle states in the vortex core, of extension O(^). These bound states 
are obtained by solving the Bogoliubov-deGennes mean-field equations for the wave 
functions in particle-hole space, in which A (a;) plays the role of the potential [78| . 
They lead to a profound alteration of the low-energy dynamical behaviour of vortices 
because, among these bound states, there exist the so-called zero modes crossing the 
zero of energy as a function of a component of the quasiparticle wave vector (the other 
core bound states excitation branches have energies at least of O(Aoo), i.e. of order 
the bulk energy gap) . These zero modes lead to a exchange of quasiparticle momentum 
between the superfluid vortex vacuum, moving with velocity dX/dt, and the quasi- 
particle heat bath, moving with velocity Vn, if the relative velocity dX /dt — t)„ is 
non- vanishing: The quasiparticles are driven by the effective (electric-like) field, stem- 
ming from this relative velocity, from the occupied negative energy levels to those of 
positive energy (energies are counted from the Fermi energy Ep), thereby transferring 
momentum to the quasiparticle bath js^ . The motion of the quasiparticles on the zero 
mode branch is a process for which the notion of spectral flow has been coined. A con- 
tribution to spectral flow is stemming only from chiral particles (cf. j72|, chapter 6). 
Due to the momentum exchange between superfluid vacuum and quasiparticle bath, 
there results an additional transverse force on the vortex, to be added to the usual 
contributions which occur in superfluids without zero mode bound states [ [77| 

The equation of motion for the vortex in the time domain is in general regions of 
the parameter space non-local. In frequency space, the equation of motion is local, 
because the convolutions of Green's functions with the vortex co-ordinates in the time 
domain become products after Fourier transformation. The equation of motion may 
then be written in a form which satisfles Galilean invariance, which implies that only 
relative velocities of line and heat bath (the lattice) as well as line and superflow are 
to occur in the equation of motion. We thus assume isotropy and the existence of 
only one charge carrier, for simplicity of representation. The force balance equation 
between Magnus and dissipative as well as transverse forces on the vortex is then 



It is assumed that the vortices are singly quantized, with a constant circulation vector, 
the factor of 1/2 in the Magnus force taking account of the paired nature of the super- 
fluid. For the point vortex we will deal with, X' — dX /da = ie^, the upper/lower 
sign valid for a positive/negative circulation vortex. The first line contains the conven- 
tional Magnus force (where ps ~ Ps{T) is the superfluid density), and the second line 
the dissipative and reactive mutual friction forces from momentum exchange between 
vortex line and quasiparticle bath. We further assume in what follows for the sake of 
simplicity that the normal component is clamped, i.e. t)„ = (this, of course, destroys 
the property of Galilean invariance afforded by (^5|)). Then, a convenient writing for 



expressed by ||7|, 1^, | 



8: 




(65) 
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the above equation of motion is 

'^Vsiuj) X X' = D{oj)X{uj) + (^D\oj) ~ ^hps^ X' X X{lj) , (66) 

that is, there is a coefficient D{uj) of the force hnear in the velocity of the vortex, the 
damping force, and a coefficient D'{oj) — hps/2 of a force perpendicular to the vortex 
velocity, which represents the Magnus force part proportional to vortex velocity (the 
Hall term). The driving term {l/2)hpsVs x X' (the superfluid current part of the 
Magnus force) is proportional to the circulation vector and the superflow current 
density = PsVs- In the case of superconductors, the driving force is equal to a 
Lorentz force exerted by the flux line on the electrons, and is given as x = 
Je\ X $0, with |$o| = h/2e equal to the flux quantum. The electrical superfluid 
current of electrons Joi = —epsVs- 

The mutual friction coefficients from equations (^5|) respectively ( |66|) are deter- 
mined by the momentum (and energy) exchange of the quasiparticle bath with the 
superfluid vacuum containing the topological defect structure vortex. In the quasiclas- 
sical limit, their value is governed by the quasiparticle kinetic equation |8(l, B£ 



4.2.1 The s-wave case 

We now consider the simple case of a singly quantized vortex in a two-dimensio- 
nal isotropic superfluid, for which the quasiparticle momentum has components k ~ 
{kr, k^), and the zero mode branch 

E{k^) = -ujok4, . (67) 

The level spacing ujq is of order h/(rn^'^) ^ A'^/Ep, where m is the (effective) fermion 
mass. The energy levels correspond to those of electrons on anomalous Landau levels, 
which are linear in momentum |8C|] . 

On the approximation level of a single relaxation-time in the quasiparticle kinetic 
equation, the coefficients are then given by |8^p| 

ihCoujo f 1 



UJ — UJo + IT ^ UJ -\- LiJt^ + IT ^ 

1 + {{loot)-^ - iu}/ujof ' 



for the longitudinal coefficient in equation (|62) and, for the transverse coefficient, 

. 1, hCQUjQ f I 1 

D [u) - -hp. 



2 '° 4 \uj — UjQ + IT ^ UJ + UJt) + IT ^ 

-hCo/2 



— ■ (69) 

1 + ((wqt)-! - lU)/uJo) 



''We remark that in the literature the setting h = 1 is frequently taken, which means that Planck's 
quantum of action h = 2n. 
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In these equations, the parameter density Cq = kp /2Tr (in three dimensions: Cq — 
kp/Sn'^) is the normal state density.[| The parameter r is a constant relaxation time 
in the collision term of the kinetic equation. The relations above hold provided that 
the conditions loq < = ksT/h, t ^ and hw < fesTc, T <^ Tc are met. 

The first condition means the thermal population of the levels represented by (|67|), 
such that the sums over the bound state levels in the general expressions for D, D' 
[ p2[ can be converted into integrals. This condition is not very stringent, and allows 
for a use of the formulas above for actually quite low temperatures, because loq ~ 
rc(Aoo(r)/Aoo(r = 0))(Aoo(r)/£;F) < r^. The second condition implies that the 
broadening of levels by scattering is much less than the temperature. Finally, the 
condition hio <C ksTc restricts the energy equivalent of the vortex motion frequency, 
fko^ to be much less than the bulk gap ~ ksTc. 

It is important to point out that the validity of ( |68|) and (|6^) relies on the transverse 
and longitudinal coefficients in the low energy limit being determined by the core level 
spacing ujq and the scattering frequency r^^ alone. The relations for the Hall and 
longitudinal coefficients are then given irrespective of an electric charge of the particles 
carrying the supcrfluid current. 

Two limits of the equations of motion are particularly well known. The first is 
provided by ujot 3> 1, lo ^ ujq and corresponds to massive vortex motion under 
influence of the superfluid Magnus force {X'^{lu) = —iojX^{Lu)): 

M,X - X' x[x- v)j . (70) 

At the low temperatures considered, we neglected the small contribution of (/i/2)(Co — 
Ps)X' xdX /dt ~ {h/2)pnX' xdX /dt, where p„ is the normal density in the superfluid 
state. This last expression represents the lordanskii force (remember that we fixed 
Vn = 0), which is also present in Bose superfluids and not related to spectral flow Jz?] . 

The Magnus force dominates in the above limit of 'slow' vortex motion with cu <ti ujq 
the massive term (the Hall term on the right-hand side is larger by a factor ofLUo/uj than 
the inertial term). Still, the vortex core dynamical mass Mc — hCo/{2LUo) ^ mCoTr^^ 
is much larger than the mass arising from the compressible superflow outside the 
vortex core. Their ratio is of order Mq/Mc ~ (d/^)^, where the quantity d signifies 
the interparticle spacing c? ^ ^, so that Mq may be neglected in the equation of motion 

Dissipative motion prevails in the limit loqt <ti 1, lot <^ 1: 

^Psv, X X' - '^CoMX + ^Co(coot)2x X X' . (71) 

The vortex motion is overdamped, with friction coefficient given by the expression 
r] = {l/2)hCQU)oT.^ It is observed that time inversion invariance is spoilt by the fact 

' The parameter Cq is in general a measure of the density at the location of gap nodes; it is zero 
if no gap nodes are present. In our case, the gap nodes are at the position of the vortex line, other 
possible occurences of nodes are in the bulk of p- or d-wave superfluids. 

j We do not show the hydrodynamic and transverse mass terms on the right-hand side of (fnl), 
which are small compared to the friction term. The transverse mass |89| , relating vortex velocity and 
momentum in different directions, can be defined dynamically from the equations of motion (pQ), like 
Mo- 
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that the first term on the right-hand side of the equation above does not have the 
factor X' , which changes sign under time inversion. The dissipation is in the given 
hmit of ohmic nature with a longitudinal conductivity independent of the driving 
frequency w. The Hall part of the Magnus force is suppressed by luot relative to the 
friction term, which is the result of the spectral flow phenomenon discussed above. 
Spectral flow is made possible because the minigap loq between zero mode levels is 
broadened by a large quasiparticle collision frequency 3> ujq. On the other hand, 
in the case of ([70|), spectral flow is impeded by the presence of the minigap (save for 
tunnelling events between the levels ||7^), and the Magnus force obtains. It is to be 
mentioned that in a superconductor the energy levels above the gap are also quantized 
into Landau levels with the inter level distance hujc, where lOq is the usual magnetic 
cyclotron frequency. For sufficiently small magnetic field, however, the spectrum is 
quasi-continous and no change of the above results applies . 



4.2.2 The d-wave case 

For p- and d-wave superconductors, the gap has nodes not only at the location of the 
vortex itself, but also in the bulk superfluid. The d-wave case distinguishes itself by 
the fact that these nodes occur on lines in momentum space, whereas in the p-wave 



case they only occur at points \ 72 . The lines of gap nodes in d-wave superconductors 
lead to a profound alteration of the effective vortex dynamical equations ||8^, |8j, |8^ . 

To explain the essential features of these changes, we observe, first of all, that 
the minigap ojq = ciJo(a) is a function of the angle a, which indicates the position of 
the gap nodes in momentum space, where the gap modulus A — Aosin(2a). The 
average minigap flo = < 0^0(0:) > is of the same order A'^/Ep as the constant luq in 
the s-wave case, ilp = 0{ujo). There exists, however, an additional energy scale, the 
true quantum-mechanical inter level distance in the vortex core Eq, which is defined 
according to a Bohr-Sommerfeld type of quantization prescription for the canonically 
conjugate quantum variables kcf,{a) and a. In a a-dependent version of (|67|), k^{a) — 
—E/ujQ{a), the Bohr-Sommerfeld quantization reads f k^{oi)da — h{ra -\- 7) (with 
7 = 1/2), so that the true interlevel distance is given by £'q"^ — h"^ da lo^^ {a) 
(the integral is rendered finite by the existence of a magnetic field 184]). There are 
several different regimes of vortex motion corresponding to the ratio of the relaxation 
rate and the frequency of vortex motion not only to the (average) minigap 0.^, 
as in the s-wave case discussed above, but also to Eq = h^/rt^uTc <C ?ii^o-0 

There is a parameter region which yields a comparatively simple and unique result. 
This is the case of ^qt <C 1, which implies i?oT ^ ft- Assuming we have in addition 
LOT <C 1, we are led back to an equation of motion of the typ e ([7l| ). On the other 
hand, if fipT ^ 1 (and ujt <C 1), there is a 'universal region' |8^, |84[], which is realized 
if Eqt <C h: The dissipative and Hall coefficients are independent of the relaxation 
time for very low temperatures and small magnetic fields and have the same order 
of magnitude. Finally, if U,qt 3> 1 and Eqt ^ /i, the equation of motion is of the 
type (|70|), with dominating Magnus force. This is a regime which will presumably 



''There is a further complication, which we do not take into account here because of the low 
temperatures T <^Tc considered. If T < Tc, the relevant energy scale approaches the cyclotron level 
spacing hwc- Furthermore, the results have in general to be written in a form taking into account 
both particle- and holelike Fermi surface parts, see |p4|. 
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not be realizable with respect to practially achieved sample purities, whereas the 
universal regime should be observable. It appears useful at this point to insert a short 
treatise on the terminology in the literature. Superconductors are classified as being 
'superclean', 'clean', 'moderately clean' and 'dirty', according to the value of the ratio 
//^, where I = vpT is the quasiparticle mean free path. The superclean limit loqt ^ 1 
(or ^Iqt ^ 1) corresponds to I 3> £,{Ep / Aao), representing a much stricter condition 
than its clean counterpart I 3> ^. The moderately clean and dirty limits correspond 
to Z ^ ^ and I ^ f , respectively. We will agree to call a superconductor 'moderately 
clean' if it has ujqt ^ 1 (but still simultaneously I ^ £.)■ Because of the above discussed 
d-wave pecularities, there is yet another notion of 'extremely clean', which ought to 
be introduced. This corresponds to the extreme limit Qqt ^ 1 and Eqt 3> fi, or 

The discussion of s-wave vortex motion to follow is thus valid for the d-wave case 
in a straightforward sense only if we are in the moderately clean region Qqt <C 1, 
Eqt <C 1, and have additionally ujt ^ 1, with the local form of the equation of motion 
in (0). 

4.2.3 Nonlocal motion in the time domain 

In order to obtain the real time motion of the vortex, we rewrite the general equation 
of motion in its convoluted form. To this end, we use that X^{u}) = —iujX'^{uj) and 
define K'^{uj) = iujD{uj), K"{u}) = iuj{D'{uj) - hp,/2). Then, by multiplying ^ 
with (l/27r) J exp(— loit), we have 



hps 
2 



-VsXX' = - J dt'K^it - t')X{t') - ^' X y dt'K"{t - t')X{t') , 
where the kernels are, by use of (pSh and (p9), the Fourier transforms 



(72) 



+ 00 



duj exp[— ici;(t — t')] 



1 + {{ujot) i-iw/wo) 

K"{t ~t') = — du; eM-Mt - t')] —— ^ ^ . . 

271- J^oo 1 + {(ujot)-^ ~ luj/uJo) 

The integrals can be evaluated via complex contour integration. The poles zi = 
loq — iT~^, Z2 = ~z\ are lying both in the lower half plane, so that the kernels are 
causal: 



K^(t-t') 
K"{t-t') 



-1 \ hCouJo 1 



I J 4 2tt 
I \ hCoUjQ 



dze-''^'~'"> 



± 



Z — Zl Z + Zl 

6{t-t') [Res(zi) + Res(-0i)] . 



1 / 4 

Calculating Res(zi) + Res(— zi) , they are given by (At = t — t'): 



K°{At) = 6{At)^ujlci^> 
K"{At) = e{At)^ujle^^ 



' At' 


( 


T 




' At' 


( 


T 





sin(ijJoAt) H cos(a;oAt) 



cos(tJoAt) sin(cjoAt) ) . 
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The equations of motion ( [70| ) and (Q), which are local in time, result from the non- 
local equations (^2|) only in case that the frequency of vortex motion is sufficiently 
smaller than the bigger of the two frequencies u>o,^t = t~^, i.e. such that to <C 
(wg + uj^Y^^ and either ujqt > 1, like in @, or ujqt <C 1, hke in (|l]), holds. 

The action corresponding to ( [72| ) is obtained by integrating with respect to the 
vortex position X{t). We arrive at 

S[X{t)] = dti^[X{t)] + J dtj dt'K^it - t')X{t) ■ X{t') 

-X' ■ J dt J dt'K"{t - t') X{t) X X{t') . (73) 

The gradient of the stream function -0 describes the superflow perpendicular to the 
vortex, Vij} = Vg x X' . 

Making use of the phase angle $ — arctan[(wo''')~^] and standard addition theo- 
rems, we can cast the action into the suggestive form 

S[X{t)] j dtij[X{t)] + ^uJ(,^ojI+ojI j dt j dt'eM-^r^t] 

X {X{t) ■ X{t') sin(tJoAi + $) - X' • (X(<) x X{t')) cos(tJoAt + $)} . 

The factors multiplying the dot- and cross-products of X{t) and X{t') in this nonlocal 
Lagrangian are for any value of loot and thus $ just 37r/2 out of phase. In case that 
ujqt ^ 1, $ ~ 7r/2, the first term with the dot-product dominates, whereas if ojqt 3> 1, 
$ ~ 0, the second one involving the cross-product does. 



4.3 Euclidean vortex motion 

For a description of tunnelling motion, we have to use the Euclidean action in the 
interval [—ft/3/2, h(3/2). Performing the Wick rotation in ( |73| ) through the replacement 
t —ite, gives the Euclidean action Se[X{te)] = —iS[X{t —ite)]- 

rhf3/2 

Se[X{te)] = -{hps/2) / dt,^[Xite)] 

J-hp/2 

+ dt, dt', [-K^'it, - t',)X{te) ■ X{t'J 

J-hp/2 J-hfi/2 
+ X^(<,-OX'.{X(t,)xX(0}] (74) 

where, under the condition that the real frequency uj <^ LUo,LUr, the kernels are ap- 
proximately given as 

f K^it,~t',) \ hCpUQ U ( Ur \ 

^ hCp I UJQ /' UJr \ 

and include nonlocality in lowest order. The dissipative kernel is of the Caldeira- 
Leggett-form ([p^- flO^ ) for ohmic dissipation, ~ (?7/7r)(Ate)~^, with the friction 
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No 


No 
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/ (ohmic) 




1 




i 


No 


No 


/ 


/ (ohmic) 


^3(9) [any e] 



Table 1 Comparison of different approaches to quantum tunnelling with respect to domi- 
nant contributions employed in the calculation of the Euclidean tunnelling action (shown by 
an entry with a checkmark, /). The Magnus column indicates if a linear coupling of the vor- 
tex velocity to the superfluid background, generating the vortex velocity part of the Magnus 
force, has been used. In the potential Vg{q), e <^ 1 means that the case of a near critical 
potential has been treated. In Ref. [p3|, this was done by an analytical method, in by a 
numerical one. In the case of entry [ [105| ] , both limitis of dominant Magnus and dissipation 
contributions have been considered, but not the intermediate case of equal strengths. This 
intermediate case has been examined numerically in |103|, far off and near criticality. 



coefficient 77 = (l/2)hCo{LUoU!r) / {loq The dissipation due to bound states is at 

a maximum if loot — 0(1) and vanishes in the limits ujot — > 00 and ujqt 0. The 
nonlocality of the Hall term in the action is of the same importance as that of the 
friction term if loq Lo-r- 



4.4 Different regions of parameter space 



It is rather obvious that the motion of a vortex for arbitrary competing contributions 
in the action ( [7^ can be quite complicated. In principle, the following contributions 
in the action are conceivable. In addition to the terms appearing in (fz^), there can be 
contributions arising from the self-interaction, like in (^, that is, the hydrodynamic 
mass and elasticity terms (the self-energy is usually absorbed into the potential) . The 
hydrodynamic mass was argued to be in general completely negligible as compared 
to the dynamic core mass in paired Fermi superfluids. The elasticity arises from 
the generalization of a 2d or rectilinear vortex to one of arbitrary shape and the 
additional self-energy this creates. Hence, in order to actually describe an imaginary 
time motion and thus evaluate the tunnelling process probability, only certain classes 
of metastability problems have been investigated. 

The (stream function) potential is, for reasons of (analytical) solvability, frequently 
represented as a quadratic plus cubic potential for a one-dimensional generalized co- 
ordinate q. This potential is conventionally parameterized with two or three quantities, 
a height h and width w of the potential, and possibly with the additional parameter 
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of closeness to criticality e.W The potential thus has the general representation 



\ 2 „ / \ 3 



(76) 



where e = — Vgj Vct measures the closeness to a critical external velocity Vcb, for 
which the barrier vanishes. A measure of the typical curvature radius of the potential 
is go- The zeros of the potential are aX q — and q = (3/2)eqo, so that the width of 
the potential may be defined to be w = {3/2)eqo. The maximum is at gmax = ego and 
its height equals h = Vge^. The different cases and approximations investigated in 
a selection of recently published papers on quantum tunnelling are brought together 
for comparison in Table |l|. For the detailed results and methods used, the reader is 
referred to the cited works. 

In relation to what we have found in the preceding section, wc can make the 
following observations. The core level spacing is of order ojq ~ fi,/(m^^), and the 
Magnus force dominates over the mass term in (|7^) if the vortex motion frequency 
Lu <^ ujQ. The collective co-ordinate approach implies that the curvature radius of the 
potential go S> ^, because the motion of a massless vortex in the potential is of typical 
frequency uj ^ h/{mqQ), the 'cyclotron' frequency associated with go- The condition 
for the collective co-ordinate approach, go ^ is thus equivalent to the dominance 
of the Magnus contribution over the core mass term in the superclean limit ujqt 3> 1. 
The fundamental hydrodynamic analysis of the last section, treating the Magnus force 
as dominant, therefore remains valid for tunnelling in the case of a superclean s-wave 
fermionic superfluid. The case of a c?-wave superfluid is, as already argued, more 
intricate. The dominance of the Magnus force only obtains in the extremely clean 
limit, whereas in a superclean limit, which has Eqt <C H, the dissipative and Hall force 
components are of comparable magnitude. This necessitates a complete treatment 
of the tunnelling phenomenon in two dimensions for this 'universal' [ |83t parameter 
regime, even for a point vortex, because, in the plane, one has to solve two coupled 
differential equations of motion. We will further discuss this case and its possible 
occurence in high-Tc superconductors below. In the extremely clean limit respectively 
for large enough cyclotron level spacing, and if the temperature goes to zero, i.e. is less 
than any of the energy scales associated with the average and true minigap flo —<ujn> 
and Eq, the superfluid Magnus force is the only remaining nondissipative force on the 
vortex, and the analysis of tunnelling in the last section is valid. 



4.5 Quantum tunnelling in high-Tc superconductors 

In conventional superconductors, the dissipative component in the vortex equations of 
motion usually is very significant, so that quantum tunnelling is largely suppressed; 
the temperature region above zero, in which temperature independent quantum tun- 
nelling is of importance, is exceedingly small |10C]. There is, however, the intensely 
investigated class of high-Tc superconductors, which can very well be in a clean or 
even superclean limit. In this latter case the Hall angle can approach the value tt/2 

'The potentia]_(r6|) corresponds to the so-called 'tilted washboard potential near criticality' if 
e < 1, see, e.g., [^OjT^ 
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(the vortex then moving more with the local superflow, rather than perpendicular to 
it) |87| . In addition, and even more important, in contrast to conventional supercon- 
ductors these materials can exhibit quite small coherence lengths in the order of the 
lattice spacing, and a ratio T* /Tc ~ 0(1 • • • 10) (cf. the energy barrier considerations 
relevant for vortex tunnelling in the introductory considerations. These facts lead to 
the possibility that in some of these superconductors, quantum tunnelling of flux lines 
might be observable at low temperatures. The ratio of the crossover temperature Tq 
from thermally activated to quantum behaviour for the flux line depinning to the crit- 



ical temperature Tc was measured in very different materials ([106|-|111|). It is found, 
depending on the material, that (To/Tc)high-Tc ~ 0.03 • • -0.09, where Tq is defined in 
these measurements to be the temperature at which fiux motion deviates from the one 
expected for purely thermal activation.^ We compare this with He II, where Tq ^ 150 
mK and ~ 2.2 K, so that (To/rc)Hoii ~ 0.07. Considering the fact that the ratio 
of the critical temperatures for these high-Tj, superconductors and helium II can be 
up to a factor of 50, the values of Tq/T'c are comparatively close, differing at most by 
a factor of two. This suggests a common physical origin of the deviation from ther- 
mal activation behaviour in helium II and high- Tc superconductors, be it quantum 
tunnelling or some other, classical, flow instability mechanism. This is only natural 
from the point of view that both systems represent, on a fundamental level, strongly 
coupled superfluids. 

There are high- Tc superconductors presumably belonging to the class of c?-wave su- 
perconductors [J74|, |7q| , or some variety of this pairing symmetry with small deviations 



from pure d-wave. The considerations of section 4.2 for the d-wave case then apply, 
provided we assume a quasiclassical, low energy treatment of vortex motion in linear 
response to be valid, at least qualitatively, in these superconductors and on the scales 
of tunnelling. 

The c?-wave superconductor, for practically achieved maximal sample purities and 
sufficiently small magnetic fields, will not be in an 'extremely clean' region (Magnus 
force dominating), but rather in the superclean 'universal' region (for the exact con- 
ditions on H and T, see [|4|). This implies that there is a magnetic field region in 
which, as already explained, the vortex tunnelling motion is not dominated either by 
the dissipative or the Magnus force (the Hall term), even for very low temperatures, 
but is governed by both forces. There are indications that in clean high- Tc materials 
an intermediate regime between purely dissipative and Hall tunnelling may indeed be 
realized ]109| , [110| . The intermediate case thus clearly needs further investigation, 
because the (measurable) temperature dependence of the Euclidean action can be dif- 
ferent from that expected for purely dissipative or Hall motion, and the transition from 
quantum to thermally activated depinning of the vortices be of first or second order 



[ pO| 111|. A first step in this direction has been made in |103|, where the problem 



of quantum tunnelling was investigated with the static (low frequency) versions of the 
formulas for the s-wave case, (^8|) and (^9|). It is found that the tunnelling rate dis- 
plays a minimum for ojqt ^ 0{\). This is in accordance with the fact that dissipation 
due to spectral fiow is at a maximum for values of the parameter uj^t which are of 
order unity. Whereas the nonlocality of the ohmic dissipation was considered there a 



In Ref. [115 1, however, quite large ratios up to (Tb/Tt.)high-Tc ~ 0.22 have been reported 
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la Caldeira-Leggett, that of the HaU term, however, was not taken into account. It 
is apparent from (75) that this is not justifiable in the intermediate case of interest, 
in which luq ~ Ur — t^^ ■ Within the formalism we presented, the Hall term can 
be treated locally in the action only in the limits loqt oo and uot 0. Apart 
from this objection, the vortex dynamical behaviour in d-wave superconductors is in 
general more complicated than in conventional s-wave superconductors, as we have 
already pointed out. For example, under certain conditions there can be resonances 
in the vortex response, if vortex tunnelling frequencies are near (2k + VjE^/h, with k 
an integer number [ p4| . A complete treatment of the tunnelling problem for a d-wave 
system in the intermediate range then necessitates an incorporation of nonlocality in 
the longitudinal and transverse vortex response, as well as possible resonances with 
collective modes induced by the moving vortex. 



5 Concluding remarks 

The present work has treated the consequences and limitations of a large scale descrip- 
tion of the motion and tunnelling generation of quantized vortices in dense superfluids. 

We may summarize the salient assumptions, pertaining to this treatment, in a 
compact way as follows. The vortex object cannot be described in its genesis^ because 
we do not know precisely in which way a vortex should be represented on curvature 
scales of order the coherence length. The intrinsic nucleation process, in all its quantum 
many-body subtleties, happens on these scales. Applying the formalism used in this 
treatise, we thus have to assert that the vortex somehow comes into a topologically 
ensured existence. We can, then, assign a collective co-ordinate to the singular center 
of topological stability. Furthermore, if we wish to describe the vortex as a string (in 
three dimensions), or point object (in two dimensions), which is tunnelling through a 
potential barrier, we have to adopt the point of view that we are allowed to quantize 
vortex position and momentum, in a canonical manner. These two assumptions and 
their validity lie at the heart of our treatment of the problem of vortex tunnelling in 
a superfluid at the absolute zero of temperature. There is, then, no theory of vortex 
nucleation in a dense, real life condensed matter system, because we drastically reduce 
the number of (quantum) degrees of freedom actually relevant for nucleation. There 
is, though, a consistent theory of vortex quantum tunnelling in the large scale domain, 
which we represented here in its formal requirements and geometric implications. 

The necessity of employing the long wavelength limit affects the contributions of 
different origin in the tunnelling exponent. The volume contribution, associated with 
the incompressible superfluid, always dominates over the area contribution. This lat- 
ter contribution is associated with the vortex effective mass, and is thus depending 
on the detailed dynamical behaviour of the vortex on the tunnelling path. The dom- 
inant volume contribution, in contrast, depends on the shape of the vortex path in 
configuration space only. 

At absolute zero, the Galilean invariancc of the bulk superfluid is required to be 
broken for tunnelling to be energetically possible. This breaking of invariance can be 
attained by considering vortex motion in the presence of obstacles. The geometrical 
implications imprinted by a specific obstacle chosen play an import part in our analysis 
and give a central result. Namely, if the vortex moves near the boundary, trapped by 
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the pinning potential generated by the obstacle, only those paths are allowed in which 
the vortex center remains at least within a distance of order the coherence length from 
the obstacle surface. If the flow obstacle, then, has curvature perpendicular to the 
flow passing at infinite distance over the obstacle much larger than parallel to this 
flow, we have shown that the (constant) energy of the tunnelling vortex cannot be 
less than a given minimal energy. This minimal energy is needed by the vortex to 
be completely describable by the collective co-ordinate with which we have equipped 
it, because else it would come within a distance less than the coherence length to 
the obstacle boundary. This result is generalizable to the case of relativistic vortices 
in spacetime, where the necessary breaking of (local) Lorentz invariancc for timelike 
currents will lead to the same kind of prediction. 

The theory we have developed can claim to make exact predictions on observable 
tunnelling probabilities in the semiclassical limit, as long as the collective co-ordinate 
approach makes (geometrical) sense and the Magnus force contribution is dominating 
the tunnelling action. We have seen, by considering equation (^3|), that the tunnelling 
scales in He II, using the available data, will be of order nanometers. Correspondingly, 
the number of particles participating in the tunnelling event, i.e. those contained in 
the volume determining the tunnelling action, is comparatively small. It is of order 
given physically realistic estimates of the prefactor and the tunnelling 
rates to be expected. Hence, the applicability of our theory, for the actual physical 
conditions encountered in He II, is restricted, in the sense that it can only give lower 
bounds for tunnelling rates, valid on large enough scales. It is, first of all, not entirely 
obvious that the dominance of the Magnus force still holds on the scales relevant for 
tunnelling. Second, we have no really direct means to compare the tunnelling rates 
observed by varying experimental conditions with the predictions of the theory. This, 
of course, stems from the very nature of the theory as a geometric theory. The pre- 
dictions it actually makes concern primarily the variation of tunnelling rates with 
the geometrical parameters of flow obstacles. Thus the only conceivable possibility 
of checking the validity of the theory is the observation of a variation of tunnelling 
rates with micro-orifice surface roughness; such an experiment, in a reproducible fash- 
ion, has not been carried out yet. In order to make further progress in relating the 
experimental findings to a suitable theory, it appears from these arguments that one 
is required to go beyond the Magnus force dominance in the tunnelling action and 
consider the modifications of this dominance in the mesoscopic scale domain, in par- 
ticular by the interaction of the vortex with the elementary excitation spectrum of 
the superfluid. This interaction should, on these scales, change the relevant effective 
forces acting on the vortex, leading to measurable effects on the tunnelling probabil- 
ity. The same considerations essentially apply to the possible observation of tunnelling 
of flux quanta {i.e., of the magnetization) in high- Tc superconductors, however with 
a considerable amount of complications, caused by different parameter domains, as 
mentioned in the third section. In addition, in the general case, the requirement that 
proper electromagnetic gauge invariancc is to be satisfied should play an important 
role. 

Future research, along the directions we have been alluding to above, is needed 
to shed further light on the intrinsic process of the genesis of quantized vortices in 
superfluids. 
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